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INTRODUCTION 


Telecommunication  is  vital  to  FAA.  The  reason  is 
more  than  just  a  means  of  communication,  but  it  is  because 
one  of  the  communication  terminals  is  usually  airborne  and 
unaccessible  by  cables.  The  growth  of  airplane  users  cre¬ 
ates  the  same  growing  demand  for  telecommunications.  All 
these  point  to  a  growing  need  for  better  spectrum  utilisa¬ 
tion  within  the  telecommunication  spectrum. 

A  control  tower  or  flight  service  station  may  commu¬ 
nicate  co  more  than  one  plane  at  a  time.  Therefore,  in 
existing  FAA  facilities,  pu! title  channels  are  provided  at 
each  location.  But  when  frequency  assignments  are  too 
close,  cross  modulation  and  inter-modulation  can  occur. 

On  one  hand,  demand  presents  a  need  to  squeeze  in  more 
charjnels  into  the  existing  spectrum.  But  on  the  other 
hand  by  having  close  frequency  assignments  it  results  in 
undesirable  cross-modulations  and  inter-modulations.  Reso¬ 
lutions  of  these  problems  require  consideration  of  the 
following: 

1.  Selectivity  of  receivers. 

2.  Sensitivity  of  receivers. 

3.  Channel  bandwidth  of  transmitters  and  power  out¬ 
put  of  transmitters . 

4.  Sharp  cutoff  bandpass  filters. 

5.  Using  RF  cancellation  techniques. 

6.  Coding  format. 

7.  Increase  antenna  isolation. 

8.  Other  techniques:  Time  sharing,  etc. 

This  report  addresses  itself  to  item  number  7 ,  im¬ 
proving  antenna  isolations. 


SOX?  FAA  COMMUNICATION  ANTENNA  SITUATIONS 

2.1  Physical  Situations 

Communication  antennas  are  sometimes  clustered 
around  a  tower  top  or  a  roof  top.  For  example,  on 
top  of  a  control  tower,  communication  antennas  share 
a  space  approximately  15  feet  by  15  feet  square. 

Thus  in  a  typical  situation,  antennas  are  mounted 
less  than  6  feet  apart. 

2.2  Existing  Communication  Antenna  Problems 

In  sites  visited,  according  to  operators,  equip¬ 
ment  performances  were  quite  satisfactory.  Occasion¬ 
ally  there,  were  some  loss  of  signal  and  fading. 

Cross  talk  is  ever  present  in  the  communication  sys¬ 
tem.  A  remedy  existing,  now  to  solve  the  problem 
rests  with  the  training  of  the  operators .  Vhev  are 
experienced  enough  to  separate  the  voice  of  the  oth¬ 
er  operator  from  the  noise  background.  The  adapta¬ 
bility  of  human  being  minimizes  the  problem.  However 
the  problem  is  still  there,  and  some  of  these  prob¬ 
lems  can  be  solved  by  better  equipment.  Loss  of  sig¬ 
nal  and  fading  can  be  due  to  antenna  problems.  They 
could  be  due  to  poor  antenna  placement,  such  that  oth 
er  antennas  or  structures  (like  lightning  rods)  inter 
feres  with  the  antenna’s  radiation  characteristics. 

It  is  seen  that  many  of  these  problems  can  be  solved 
or  minimized  by  improving  antenna  isolations . 


SOLUTIONS  TO  THE  ANTENNA  ISOLATION  PROBLEM 

When  one  looks  at  the  problem  of  increasing  isola¬ 
tion,  and  if  one  takes  a  systems  approach  one  must  first 
list  the  antenna  parameters  that  will  affect  isolation. 

The  maior  parameters  involved  here  are: 

1.  Amplitude, 

2 .  Phase ,  and 

3.  '  Polarization. 

These  are  the  electrical  parameters  to  be  consid¬ 
ered  and  of  course  there  are  mechanical  parameters  like: 

1.  Size, 

2.  Rigidity, 

3 .  Weatherability  , 

Reliability.  And  lastly  but  not  least, 

5.  Compatibility  with  the  rest  rf  the  communication 
system. 

Consider  first  the  three  major  electrical  parameters, 
they  are  taken  up  one  at  a  time  in  the  following  sections : 

3 . 1  Amplitude  Consideration 

In  this  topic  one  takes  into  consideration  an¬ 
tenna  spacing  and  pattern  shape.  For  FA/'  communica¬ 
tion  facilities,  the  desired  pattern  shape  is  gener¬ 
ally  omnidirectional  in  the  azimuth  plane.  Therefore 
isolation  by  taking  advantage  of  pattern  is  usually 
not  applicable.  The  only  other  thing  is  antenna  spa" 
ing.  To  achieve  30  -  40  db  isolation,  the  two  anten¬ 
nas  have  to  be  spaced  80  or  so  feet  apart.  This  is 
of  course  not  possible  if  one  confines  the  m irennas 
to  a  common  roof  top  or  tower. 

3.2  Phase  Considerations 

Existing  ideas  involve  the  use  of  a  feedback 
loop.  If  the  feed  back  contains  frequencies  and 
phases  that  are  different  from  the  desire  i  signal, 
these  extraneous  signals  can  be  suppressed  by  phase 
opposition.  It  is  difficult  to  determine  and  use  the 
coupling  characteristics  of  all  the  antennas  involved 
in  order  to  create  proper  phase  cancellation,  without 
deteriation  of  the  desired  signal. 


Another  method  without  using  the  feedback  1 oop 


ends  up  as  a  circular  array  system.  The  system  is 
discussed  in  section  4 . 

3.3  Polarization  Considerations 

Isolation  by  polarization  diversity  is  very  ef¬ 
fective.  But  for  air-ground  communication  in  use  by 
FAA,  the  types  of  polarizations  applicable  are  limit¬ 
ed.  They  are  limited  to: 

1.  Vertical  Polarization, 

2.  Right  Handed  Circular  Polarization, 

3.  Left  Handed  Circular  Polarization. 

Polarization  diversity  using  a  vertically  polar¬ 
ized  antenna  and  a  horizontally  polarized  antenna  can 
not  be  used  in  FAA  facilities.  Another  way  is  to  use 
circular  polarizations  of  different  rotations.  It  has 
been  shown  that  relatively  high  isolation  can  be  a- 
chieved  by  using  such  an  arrangement.  However,  this 
technique  only  provides  one  additional  isolated  antenna. 

3.4  Vertical  Stacking 

By  vertically  stacking  two  vertically  polarized  an- 
jennas ,  in  theory,  high  isolation  can  be  achieved.  But 
in  practice,  it  is  impossible  to  achieve  high  isolation 
if  existing  antennas  are  used,  because  the  feeding  cable 
of  the  antenna  on  top  will  interfere  with  the  radiation 
characteristics  of  the  lower  antenna.  Furthermore,  iso¬ 
lation  value  is  very  sensitive  to  the  alignment  of  the 
two  antennas.  However  this  does  not  exclude  the  possi¬ 
bility  of  designing  two  antennas  coaxially  mounted. 

Such  that  the  feed  cable  of  the  antenna  on  top  can  go 
straight  through  the  body  of  the  lower  antenna. 

3.5  Some  Mechanical  Considerations 

Some  desirable  features  that  may  be  incorporated 
into  a  design  are  that  the  antennas  should  be  construc¬ 
ted  as  a  module  so  that  stacking  can  be  done  to  in¬ 
crease  the  number  of  units.  Some  kind  of  weather  pro¬ 
tection  shield  should  be  provided  .in  a  form  of  radomes 
to  minimize  weather  interference,  and  also  to  hold  an¬ 
tenna  tolerances.  The  construction  of  supporting  struc¬ 
tures  should  be  designed  so  as  to  minimize  interference 
with  antennas'  electrical  characteristics.  For  example, 
tower  railings  can  create  a  ground  plane  effect.  And 
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so,  lightning  rods  can  shadow  the  antenna  radiation, 
is  impossible  to  design  an  effective  antenna  system 
without  taking  into  consideration  of  supporting  struc¬ 
tures. 


4.  CIRCULAR  ARRAY  CONCEPT 


It  is  a  known  fact  that  if  a  circular  array  is  fed 
with  a  progressive  phase,  it  would  create  a  null  at  the 
array  center.  If  an  antenna  is  placed  at  the  null  point, 
then  perfect  isolation  can  be  achieved.  This  is  the  con¬ 
cept  of  using  circular  array  to  achieve  antenna  isolation. 
Further,  if  the  array  is  fed  by  a  Shelton-Butler  matrix 
tnen  multiple  beams  can  be  formed.  In  this  case,  each  de¬ 
sired  beam  is  doughnut  shaped  and  is  formed  by  a  different 
phase  progression. 

4.1  Advantages  of  Using  a  Circular  Array 

Because  of  the  possibility  of  simultaneous  beam 
forming,  many  isolated  channels  equi valent  to  many  i- 
solated  antennas  can  be  formed  simultaneously.  Axi- 
other  advantage  is  the  small  size,  because  a  circular 
array  forms  the  isolated  channels  within  itself,  the 
size  is  thus  limited  to  the  size  of  the  circular  array 
itself.  High  degree  of  isolation  is  achieved  through 
precision  and  not  by  size. 

The  concept  of  the  phased  circular  array  is  the 
mixture  of  amplitude  consideration  and  phase  consider¬ 
ation,  in  that  all  the  array  elements  must  be  fed  with 
equal  amplitude  and  the  phase  is  progressive.  If  the 
amplitudes  are  not  equal,  then  the  z'esultant  null  would 
not  be  a  deep  one.  Using  this  concept,  polarization 
diversity  does  not  have  to  enter  into  the  picture. 
However,  polarization  may  be  used  to  further  the  advan¬ 
tage  gained  by  using  the  array. 

The  array  analyzed  in  this  report  employs  verti¬ 
cally  polarized  elements .  The  concept  can  be  extend¬ 
ed  to  circularly  polarized  antennas , 

In  order  to  achieve  perfect  isolation  the  element 
to  be  isolated  from  the  array  must  be  placed  perfectly 
at  the  null  point.  But  in  fact  that  is  not  possible, 
therefore  a  sizable  portion  of  this  report  is  devoted 
to  finding  the  sensitivity  of  position  error  on  the 
magnitude  of  isolation.  In  the  next  section  the  theo¬ 
retical  formulation  of  the  sensitivity  study  will  be 
presented  and  its*  results  discussed. 
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5 .  THEORETICAL  STUDIES 


The  basic  technique  employed  in  this  research  is  the 
moments  method.  The  formula  used  is  the  Haller* s  Integral 
Equation.  The  array  configuration  analyzed  is  shown  in 
Figure  1.  A  computer  program  is  written  to  analyze  such 
an  array,  in  the  analysis  it  is  assumed  that  there  are 
four  elements  equally  spaced  around  a  circle  and  a  fifth 
element  placed  near  the  center  of  the  circle.  The  loca¬ 
tion  of  the  center  element  is  varied  to  study  its  effect 
of  radiation  patterns,  current  distributions,  and  isolation. 
The  computer  program  itself  however  is  quite  general.  It 
can  handle  array  elements  placed  in  any  location  and  they 
do  not  have  to  be  placed  on  a  circle . 

The  center  element  is  loaded  with  a  fixed  resistance. 

In  the  program  the  value  of  the  resistance  can  be  varied, 
but  for  the  purpose  of  this  study  and  for  uniformity  of 
results  in  order  to  facilitate  comparisons,  the  resistance 
value  is  fixed  at  100  ohms. 

The  detailed  discussion  of  the  Moments  Method,  the 
formulation  of  Ha lien's  Integral  can  be  found  in  Harrington 
[1]  and  Kraus  [2].  Assumptions  employed  in  writing  the  com¬ 
puter  program  and  the  program  listing  are  included  in  Appen¬ 
dix  A, 

The  dipole  elements  are  assumed  to  be  half  wavelength 
long  center  fed  and  symmetrical .  The  feed  points  are  all 
on  the  same  plane.  The  feed  voltages  are  labeled  V^,  V<2  to 
Vg  as  shown  in  Figure  1.  A  unique  technique  called  sequence 
function  method  is  employed  while  calculating  the  input  im~ 
pecances  and  the  technique  is  described  in  Appendix  B.  The 
output  of  the  program  shown  in  Appendix  A  includes  the  gen¬ 
eralized  admittance  matrix,  which  leads  to  current  distribu¬ 
tions  and  input  impedances  of  the  antennas. 

A  separate  program  using  a  different  approach  also  cal¬ 
culates  the  current  distribution  and  impedances.  The  pro¬ 
gram  is  listed  in  Appendix  C, 

A  third  program  which  uses  the  output  of  the  second 
program  to  calculate  input  impedances  and  current  distribu¬ 
tions  along  with  a  fourth  program  which  calculates  radia¬ 
tion  patterns  and  phases,  as  well  as  plotting  th«_m,  are  also 
included  in  Appendix  C. 

A  fifth  program  which  uses  the  current  distribution 


7 


da-ra  to  calculate  isolation  data  is  appended  in  Appendix  C. 
A  sixth  program  which  scales  and  superpositions  the  phase 
sequences  to  obtain  final  field  patterns  is  also  appended 
in  Appendix  C,  Theoretical  results  have  been  obtained 
using  the  four  phase  progressions  shown  in  Table  1.  They 
are  discussed  in  following  sections.  It  is  assumed  that 
the  center  element  is  terminated  with  a  100  ohm  load. 

f«.l  Current  Distributions 

After  obtaining  the  generalized  admittance  ma¬ 
trix,  the  first  information  that  is  measurable  is  the 
current  distribution  along  the  length  of  the  dipole 
antennas.  The  case  where  the  array  diameter  is  0.3 
wavelength  and  the  di jole  cross-section  radius  0,025 
wave-length  is  considex'ed  first.  Figures  2  through  4 
show  current  distributions  for  four  different  phase 
sequences.  It  is  seen  in  Figure  2  that  the  current 
of  the  center  element  for  the  zero  degree  progression 
is  very  large.  That  is  because  all  the  dipole  ele¬ 
ments  are  in  phase  and  their  powers  reinforce  each 
other  at  the  location  of  the  center  element.  However, 
Figure  3  shows  the  current  distribution  of  the  case  of 
plus  ana  minus  90  degrees .  It  is  seen  there  that  the 
center  element  currant  distribution  is  very  small  and 
that  is  the  basis  of  high  isolation  when  the  phased 
circular  array  concept  is  used.  Figure  4  shows  the 
similar  current  distribution  when  the  phase  progres¬ 
sion  is  increased  to  180  degrees.  The  current  dis¬ 
tribution  on  the  center  element  is  again  very  small. 

Results  of  other  array  confirmations  are  shown 
in  Appendix  D. 

5.2  Short-Circuited  Admittances 

Tables  2  through  4  show  the  short-circuited  ad¬ 
mittances  of  the  above-mentioned  array.  They  are  pro¬ 
vided  as  an  intermediate  step  for  other  calculations. 
Namely,  input  impedance  and  admittance  calculations 
and  load  current  calculations . 

Results  of  other  array  configurations  are  shown 
in  Appendix  D. 

5 . 3  Input  Impedance 
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PHASE  PROGRESSIONS 


0° 

90° 

-90° 

oo 

o 

0 

1201 

1/QI 

1/DI 

1/01 

1/QI 

1/301 

12=201 

1/1801 

1201 

1/Ml 

1/=1S& 

1ZQ1 

m i: 

1Z2Z01 

12=2Zfll 

1ZM°, 

2%  POSITION  ERi 
0°  PROGRESSION 

- —REAL 

- IMAGIN 


n  POSITION  ERROR 


90°,  -90°  PROGRESSION 

- REAL 

-  IMAGINARY 

(  )  -90° 


Figure  3  :  Current  Distribution  for  the  +90°  and  -90 
Progressions.  Radius  of  the  Array  (R)  is 
of  the  Dipole  (A)  is  0.025  A>  and  Position 
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CURRENT  (in  mA) 
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Figure  4  :  Current  Distribution  for  the  180°  Phase  Progression. 
Radius  (R;  of  the  Array  is  0.3  X,  Radius  of  the 
Dipole  (A)  is  0,025  X,  and  Position  Error  is  2%. 
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A  =  0.0250 
R  =  0.3000 
X5=  0.0015 
Y5=  0.0000 

ERR~R  =0.5% 

Yll  =  0.6166lE-02+j0 
Yi3  =-0. 67795E-03-j0 
Y15  =  0.74757E-04+j0 
Y22  =  0. 61705E-02+j0 
Y24  =-0. 67797E-03-j0 
Y  3 1  =-0.  67795E-03-j0 
Y33  =  0. 61749E-02+j0 
Y35  =  0 . 72273E-04+ jO 
Y42  =-0. 67797E-03-j0 
Y44  =  0. 61705E-02+j0 
Y51  =  0.74758E-04+j0 
Y53  =  0 . 7227  6E-04+j0 
Y55  =  0 . 45548E-03- jO 

Zn  =  0 . 12410E+C3+ jG 
Z13  =-0. 15344E+02+jO 
Z15  =  0 . 79162E+01- j  0 
Z22  =  0 . 12437E+03+ jO 
Z24  =-0  1534  3E+02+j0 
Z31  =-0. 15344E+02+j0 
Z33  =  0 . 12464E+03+ jO 
Z35  =  0. 68464E+01-j0 
Z42  =-0 . 15343E+02+ jO 
Z44  =  0. 12437E+03+j0 
Z  Si  =  0 . 79162E+01- j  0 
Z53  =  0.68463E+01-j0 
Z55  =  0 . 86780E+02+j0 


.  99869E-03 
.4267  2E-03 
.  41536E-02 
.  10275E-02 
. 42664E-03 
.  42672E-03 
. 10559E-02 
.  41070E-02 
. 42664E-0  3 
. 10275E-02 
.4153SE-02 
, 41070E-02 
.  15182E-02 

. 91576E+01 
. 38775E+02 
.  53  309E+02 
, 91088E+01 
.  3877  6E+02 
, 38775E+02 
. 90552E+01 
. 52459E+02 
. 38776E+02 
. 91089E+01 
. 53309E+02 
. 52459E+02 
.  557  20E+02 


0 .90292F 
0.90292E 
0.90292E 
0.9C732E 
0 .73529E 
0.90732E 
0.S0732E 
0.90292E 
0.S0732E 
0.73529E 
0.73532E 
0.73532E 


03+j  0 
03+j0 
03+ jO 
03  + jO 
04+j0 
03  + jO 
03  +j  0 
03+ j  0 
03+i0 
04+ jO 
04+ jO 
04  + jO 


.12915E-02 
.  12915E-02 
.  12915E-02 
. 13201E-02 
. 41301E-02 
. 13201E-02 
. 13201E-02 
.12915E-02 
.13201E-02 
.  41301E-02 
.41301E-02 
.  41301E-02 


-0 .50212E+02 
-0 . 50212E+02 
-0 . 50212E+02 
-0 .49943E+02 
0  3753E+01 
-0 . 49943E+02 
-0 .49943E+02 
-0 . 50212E+02 
-0 . 4994  3E+02 
0.73753E+01 
0 . 737 53E+01 
0 .73753E+01 


-jO .  33150E+02 
-jO .  33150E+02 

-  jO  .33150E+02 

-  jO . 33  201E+02 

-  jO .  52885E+02 
-jO .  33201E+02 

-  jO .  33201E+02 

-  jO . 33150E+02 
-jO  .33201E+02 
-jO  .52885E+02 
•  j  0 . 52885E+02 
-jO . 52885E+02 


Table  2:  Short  Circuit  Admittance  (Y’s)  and  Open  Circuit 
Impedance  <7's)  of  Circular  Array  Elements  and 
the  Center  Llement.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 
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A  =  0.0250 
R  =  0.3000 
X5=  0.0060 
Y5=  0.0000 

ERROR  =2.0% 

Ylx  =  0 . 61527E-02+ jO . 91060E-03 
Yi3  =-0 . 67772E-03- jO . 42652E-03 
Y15  =  0.78434E-04+j0.42251E-.02 
Y22  =  0. J1705E-92+j0.10289E-02 
Y  24  =-0. 67805E-03-j0.42529E-03 
Y31  =-0 . 67772E-O3-j0 .42652E-G3 
Y33  =  0 . 61879E-02+j0 . 11396E-02 
Y35  =  0.68  518E-04+i0 .40388E-02 
Y42  =~0 . 67805E-03- jQ . 42529E-Q3 
Y44  =  0 . 61705E-Q2+j0 . 10 289 E- 02 
Y51  =  0 . 7843  3E~04+j0 . 42251E-02 
Y53  =  0 . 68520E-04+j0 .40388E-02 
Y55  =  0 . u5S56E-03- jO . 15211E-02 

Z11  =  0.12327E+03+j0.927  64E+01 
Z13  =-0.15330E+02+jG.38771E+02 
Z15  =  0 . 95590E+01- jO . 54543E+02 
Z22  =  0 . 12439E+03+j  0 . 91139E+01 
Z24  =-0 . 15324E+02+j0 .38781E+02 
Z31  =-0 . 15330E+02+ jO . 38771E+02 
Z33  =  0 . 1254lE+03+j0 . 88690 E+01 
Z35  =  0.52809E+01-j0.51144E+02 
Z42  =  0 .15324E+02+j0 .38781E+02 
744  =  Q . 12439E+Q3+ jO . 91139E+01 
Z51  =  0 . 95590E+01- j 0 .54543E+02 
Z53  =  0.52809E+01-j0 . 51144E+0  2 
Z55  =  0,86762E+02+j0.55632E+02 


Y12  =  0 . 89628E~03+j0 . 12488E-02 
Y14  =  0 . 89628E-03+30 . 12438 E- 02 
Y2i  =  0 . 89628E-03+j  0 . 1248  8E-02 
Y23  =  0 . 9i387E-03+ jO . 13632E-02 
Y25  =  0 . 73740E-04+ j  0 .41299E-02 
Y32  =  0 . 913 8 7E- 03+10 . 13632E-02 
Y34  =  0 . 91387E~03+j0 .13632E-02 
Y41  =  0. 89628E-03+j0.12488E-C2 
Y43  =  0. 91387E-03+j0 . 13632E-J2 
Y45  =  0 . 73740E-04+ j  0.41299E-02 
Y52  =  0 . 73740E-04+T 0 ,4129^E*02 
Y54  =  0.73740E-04  +  j0.4129C'-:; 


Z12  =-0. 506Q9E+02~i0  .33080E+02 
Z14  =-0.50609E+02-j0 . 33080E+02 
Z21  =-0.50609E+02-j0.33080E+02 
Z23  =-0 .  49  53  5F.+02- jO .33285E+02 
Z25  =  0 . 73234E+01-i0 . 52867E*02 
Z32  =-0.4953  5E+02- jO . 3328  5E+02 
Z34  =-Q ,49535E+02~j0 .33285E+02 
Z41  =- 0.50 609 E+02-j0.33080E+02 
Z43  =-0.49535E+02-j0  .33285E+02 
Z45  =  0 . 73234E+01- jO . 52865E+02 
Z52  =  0,73234E+01~j0.52865E+02 
Z54  =  0 . 73234E+01- jO . 52865E+02 


Table  3:  Short  Circuit  Admittance  (Y’s)  and  Open  Circuit 
Impedance  (Z’s)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 
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A  =  0.0250 
R  =  0.3000 
X5=  0.0300 
Ys=  0.0000 

ERROR  =  10.0% 

YU  =  0.S0766E-02+30.39141E-03 
Y13  =-0 . 671S9E-03~ jO . 42157E-Q3 
Y15  =  0 . 96820E-04+ jO . 46448E-02 
Y22  =  0. 61686E-02+j0.10635E-02 
Y24  =-.Q .  67997E-03- jO .  39072E-03 
Y31  =-0 . 67189E-03- jO . 42157E-03 
Y33  =  0 . 62516E-02+jQ . 15457E-02 
Y35  =  o . 4812lE~04+j0 . 37055E-02 
Y42  =-0 . 67997E~03~ jO .39072E-03 
V44  =  0. 61686E-02+j0. 10635E-02 
Y51  =  0 . 96819E-04+j0 . 46448E-02 
Y53  =  0 . 48122E-04+j0 . 37055E-02 
Y55  =  0.4577  6E-03- jO . 15965E-02 

Zii  =  0.11829E+03+j0.9097lE+01 
Z13  =-0. 14988E+02+j0 .38652E+02 
Z15  =  0 . 18976E+02- jO . 6G223E+02 
Z22  =  0 ,12490E+03+j0 .92330E+01 
Z-?4  =-0 . 148  22E+02+j0 . 38900E+02 
Z31  =-0 . 14988E+02+j0 . 38652E+02 
Z3'-i  =  0 .12886E+03+j0  .73442E+01 
Z35  =-0 . 2186lE+01-j0 .43429E+G2 
Z42  =-0 . 14822E+02+j0 . 38900E+02 
Z44  =  0 . 12490E+03+j0 .92330E+01 
Z51  =  0. 18976E+02-j0 .60223E+02 
Z53  =-0.21862E+01-j0.43429E+02 
Z55  =  0 .86336E+02+j0 . 53419E+02 


Yl2  =  0 . 86030E~03*jO . 10255E-02 
Ym  =  0 . 88030E-03+ jO .  10255E-02 
Y2i  =  0 . 86G30E-03+ jO .10255E-02 
Y03  =  0 . 94770E-03+ jO .15956E-02 
Y25  =  0 . 79073E-04* jO .41244E-02 
Y32  =  0 . 94770E-03+ jO .15956E-02 
Y34  -  0 . 94770E-03+ jO . 1595SE-02 
Y41  =  0 . 86030E-03+ jO .10255E-02 
Y43  =  0 . 94770E-03+ jO . 15956E-02 
Y45  a  0 . 79073E-04+ jO . 41244E-02 
Y52  =  0.79073E-04+j0. 41244E-02 
Y54  =  0 .79073E~04+i0 .41244E-02 


Z12  =-0 . 52565E+02- jO . 32794E+02 
Zi4  =-0 . 52565E+02- jO . 32794E+02 
Z2i  =~0 - 52565E+02- jO . 32794E+02 
Z23  =-0 . 47305E+02- jO . 33921E+02 
Z25  =  0 . 600  56E+01- jO . 523  55E+02 
Z32  =-0.47305E+02-j0 .33921E+02 
Z34  =-0 . 47305E+02- jO . 33921E+02 
Z41  =-0 . 52565E+02- jO .32794E+02 
Z43  =-0 . 47305E+02- jO . 33921E+02 
Z45  =  0 . 60056E+02- jO . 52355E+02 
Z52  =  0 . 60056E+01- jO . 523  55E+02 
Z54  =  0 . 60056E+02- jO . 52355E+02 


Table  4:  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 

Impedance  (Z's)of  Circular  Array  Elements  and  the 
Center  Element.  Values  Shown  Follow  the  Order: 
Real,  and  Imaginary. 
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It  is  desirable  to  knoc7  how  sensitive  the  array's 
input  impedance  is  to  the  center  element's  position. 

This  gives  the  information  on  how  the  feeding  network 
must  be  adjusted  to  achieve  good  power  transfer. 

For  the  0°  progression,  and  when  array  radius  is 
0.3  wavelength,  it  is  seen  in  Figure  5  that  the  real 
part  of  the  impedance  (resistance)  does  not  change 
with  position  error,  and  some  change  in  the  imaginary 
part  (reactance)  as  shown  in  Figure  6.  Ports  1  and  3 
are  affected,  because  the  center  element  is  displaced 
towards  element  1  and  away  from  element  3,  whereas  the 
distance  from  elements  2  and  4  are  practically  unchanged. 

When  the  array  radius  is  increased  to  0.5  wave¬ 
length,  results  are  shown  in  Figures  7  and  8.  Resist¬ 
ances  of  ports  1  and  3  are  seen  to  change  with  position 
error . 

Returning  to  the  case  where  array  radius  is  0.3 
wavelength,  but  changing  the  phase  progression  to  90°, 
it  was  found  that  all  ports  of  the  circular  array  are 
insensitive  to  position  error.  Results  are  shown  in 
Figures  9  and  10.  However,  it  is  seen  in  the  figures 
that  the  input  impedance  of  the  center  element  (port  5) 
changes  drastically  with  position  error,  and  seems  to 
approach  infinity  at  zero  position  error.  That  means 
for  a  fixed  load  at  that  port ,  impedance  mismatch  is 
greater  when  the  center  element  is  closer  to  array  cen¬ 
ter.  In  the  same  manner,  isolation  between  the  element 
and  the  array  goes  up  with  decreasing  error. 

5.4  Radiation  Pattern 

Calculations  were  made  of  the  circular  array  where 
the  center  element  was  loaded  with  100  ohms.  Results 
are  plotted  in  the  following  figures. 

Figures  11  through  13  are  for  an  array  with  0.3 
wavelength  radius  at  increasing  position  errors.  It 
is  seen  that  pattern  minima  of  the  0°,  90°,  180°  phase 
progressions  are  -1.2,  -3.2  and  -codb  respectively. 
Patterns  of  the  +90°  and  -90°  progressions  are  identi¬ 
cal.  With  increasing  error,  pattern  distortions  are 
seen  to  be  within  1  db . 

Figures  14  through  16  show  radiation  patterns  of 
the  array  when  array  radius  is  increased  to  0.5  wave- 
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RESISTANCE  | 


POSITION  ERROR 


Figure  5  :  Magnitude  of  Input  Resistance  as  a  Function,  of  Position 
Error  for  the  0°  Progression,  0.3  Wavelength  Radius. 
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POSITION  ERROR 

Figure  7  :  Magnitude  of  Input  Resistance  as  a  function  of  Position 
Error  for  the  0°  Progression,  0.5  Wavelength  Radius. 


POSITION  ERROR 


Figure  9 


:  Logarithmic  Input  Resistance  as  a  Function  of  Position 
for  The  90®  Progression  0.3  Wavelength  Radius. 


RELATIVE  POWER  (db) 


A  =  0,0250  A  R  =  0,3000  A 
X5=  0,0060  A  Y5=  0,0000  A 
ERROR  =  2,6% 


Figure  12 


$  -  PLANE  ANGLE  (Degrees) 

Theoretical  Azimuthal  Plane  Radiation  Pattern 
of  the  Circular  Array. 
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A  =  0.0250  A  R  =  0.5000  X 

x5=  0,0100  X  y5=  0.0000  A 

ERROR  =  2,0? 


RELATIVE  POWER  (db) 


A  =  0,0250  A  R  =  (J ,5000  X 
X5=  0,0500  A  Y5=  0,0000  A 
ERROR  =  10,0% 


js  -  PLANE  ANGLE  (Degrees) 

Figure  16  :  Theoretical  Azimuthal  Plant.  Radiation  Pattern 
of  the  Circular  Array. 


29 


length.  It  is  seen  that  lobes  are  formed  at  all  phase 
progressions.  As  error  increases,  patterns  are  dis¬ 
torted  in  shape.  The  array  with  0.5  wavelength  radius 
is  thus  seen  unsuitable  for  omnidirectional  communica¬ 
tions  . 

Figures  17  and  18  show  patterns  of  the  0.3  wave¬ 
length  radius  array  as  in  Figures  11  through  13  but 
the  displacement  of  the  center  element  is  in  the  di¬ 
rection  between  two  adjacent  array  elements  (0  -  45°) 
rather  than  towards  one  element  (0  =  0  ° ) .  They  are 
similar  to  those  shown  in  Figures  11  through  13. 

Other  patterns  including  those  with  a  different 
dipole  radius  are  included  in  Appendix  D. 

5.5  Isolation 

Isolation  is  defined  as  the  ratio  of  power  trans¬ 
ferred  to  the  100  ohm  load  to  the  total  power  entering 
the  array . 

Power  transferred  to  the  100  ohm  load  is  equal  to 
the  product  of  the  100  ohm  and  the  square  of  the  feed 
point  current  of  the  center  element.  Power  entering 
the  array  is  the  sum  of  powers  entering  each  array  ele¬ 
ment.  Each  array  element  is  assumed  to  have  a  unit 
voltage  across  its  feed  point.  Computer  program  for 
isolation  can  be  found  in  Appendix  A. 

Results  are  shown  in  Figures  19  and  20,  Figure  19 
shows  isolation  as  a  function  of  position  error.  It  is 
seen  that  except  for  the  0°  progression,  all  other  three 
progressions  show  increasing  isolation  with  decreasing 
error.  It  approaches  infinity  asymptotically.  The  re¬ 
sult  verifies  the  trend  of  input  impedance  discussed  in 
section  5.3. 

Figure  20  shows  the  effect  of  dipole  radius  on  iso¬ 
lation.  For  the  range  of  radii  investigated,  isolation 
is  insensitive  to  dipole  radius. 
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Figure  20  :  Isolation  as  a  Function  of  Dipole  Radius. 
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6. 


EXPERIMENTAL  STUDY 


A  laboratory  setup  was  assembled  using  scaled  models 
of  existing  FAA  antennas  and  a  power  dividing  and  phasing 
network.  This  section  is  divided  to  cover  antenna  and  net¬ 
work  design*  radiation  characteristics,  isolation,  gain, 
and  other  capabilities  of  the  array. 

6.1  Antenna  and  Network  Design 

An  FAA  coaxial  dipole  was  used  as  the  basis  of 
design.  The  FAA  dipole  was  originally  designed  for 
use  at  VHF.  Ey  keeping  the  proportion,  the  dipole 
was  scaled  down  in  size  to  L-band,  with  a  center  fre¬ 
quency  of  1.45  GHz.  The  mechanical  drawing  of  the 
scaled  antenna  is  shown  in  Figure  21.  An  exploded 
view  of  the  scaled  dipole  is  shown  in  Figure  22 .  The 
assembled  dipoles  are  shown  in  Figure  23. 

VSUR  of  dipoles  were  measured.  VSWR  plot  of  the 
original  VHF  dipole  is  shown  in  Figure  24 ,  and  the 
plot  of  scale  models  is  shown  in  Figure  25.  It  is 
seen  that  general  shapes  of  the  plots  are  similar, 
with  the  exception  that  the  VHF  dipole  has  a  much 
narrower  bandwidth  but  a  better  match  at  center  fre¬ 
quency.  This  is  because  the  VHF  antenna  tested  had  a 
matching  network  at  the  feed  point.  The  matching 
network  of  the  original  dipole  consists  of  a  single 
coil.  Because  of  that  matching  network,  the  antenna 
is  tuned  to  a  single  frequency  and  exhibits  a  nar¬ 
rower  bandwidth  characteristic  there.  Since  a  per¬ 
fect  match  is  not  required  to  prove  the  theory  of  the 
circular  array ,  the  matching  network  was  not  used  for 
the  scale  model . 

For  ease  of  controlling  phase  distribution,  a 
simple  power  divider  system  incorporated  with  varia¬ 
ble  phase  shifters  was  used.  The  circuit  diagram  is 
shown  in  Figure  26.  The  input  power  is  fed  in  from 
the  top  and  is  first  divided  into  two  equal  parts  by 
an  isolated  power  divider  then  further  divided  into 
four  parts  by  two  more  isolated  power  dividers.  At¬ 
tached  to  each  of  the  four  arms  of  the  divider  is  a 
trombone  line  stretcher.  After  the  line  stretcher, 
the  line  feeds  into  an  antenna.  Figure  27  shows  the 
picture  of  such  a  network. 

The  accuracy  of  phase  adjustment  was  within  ±5° 
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Figure  21 


The  Mechanical  Drawing  of  The  Scaled 
Coaxial  Dipole . 
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FREQUENCY  (MHz) 

Figure  24:  VSWR  of  VHF  Dipole  Type  FA-52^8. 


39 


ISOLATED  POWER 
DIVIDERS 


VARIABLE 

MECHANICAL 

PKASESHIFTERS 


TO  CIRCULAR  ARRAY 


Figure  26:  Circuit  Diagram  of  the  Feed  Network. 


mmm. 

and  the  accuracy  of  power  distribution  was  within 
±1  db. 

A  circular  array  with  a  diameter  of  0.6  wave¬ 
length  was  formed  using  above  described  elements 
and  feeding  network.  The  array  itself  was  supported 
by  rigid  dielectric  braces.  The  trombone  phase 
shifters  were  adjusted  to  give  a  90°  phase  progres¬ 
sion.  The  antenna  feed  arrangement  is  depicted  in 
Figure  28,  where  element  No.l  is  fed  with  unit  am¬ 
plitude  and  0°  phase,  element  No. 2  is  fed  with  unit 
amplitude  and  90°  phase,  element  No. 3  is  fed  with 
unit  amplitude  and  180°  phase,  and  the  last  one  is 
fed  with  270°  phase.  The  element  located  at  the 
center  is  the  test  element  which  is  supposed  to  be 
isolated  from  the  circular  array  by  virtue  of  the 
array’s  phase  progression.  The  system  was  then 
tested,  and  results  discussed  in  following  sections. 

6.2  Radiation  Characteristics 

The  array  was  tested  for  three  different  phase 
progressions .  Array  patterns  and  center  element  pat¬ 
terns  were  plotted  as  well  as  other  variations.  They 
are  discussed  below. 

6.2.1  The  90°  progression 

Radiation  patterns  of  the  90°  progression  are 
shown  in  Figure  29 .  The  array  pattern  has  a  peak  to 
null  variation  of  +3.5  db.  But  part  of  that  varia¬ 
tion  was  caused  by  network  amplitude  unbalance,  which 
shows  up  as  a  wave  in  the  pattern.  If  the  network  is 
balanced  and  the  wave  is  taken  out  from  the  pattern, 
it  will  leave  only  the  ripples ,  and  the  variation 
will  be  reduced  to  +2  db. 

Because  of  the  array  blockage,  the  pattern  of 
the  center  element  is  seen  to  have  ripples  too.  The 
pattern  is  shown  in  the  same  figure.  But  it  is  seen 
that  the  pattern  has  less  variation  than  the  array 
pattern. 

When  the  center  element  position  error  is  in¬ 
creased  from  2%  to  10%,  it  is  found  that  the  array 
pattern  stays  practically  the  same,  but  the  pattern 
of  the  center  element  has  changed  in  two  ways.  First 
of  all  the  relative  amplitude  has  decreased  in  some 
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spots ,  and  then  the  pattern  is  less  isotropic .  The 
patterns  are  shown  in  Figure  30.  The  reason  for  the 
changes  is  because  the  array  is  now  coupled  to  the 
center  element,  closer  than  before,  due  to  an  increased 
position  error  of  the  center  element. 

In  a  reverse  manner ,  the  pattern  of  the  center 
element  will  be  more  isotropic  and  has  a  higher  gain 
if  coupling  is  less.  This  will  be  shown  to  be  true  in 
section  6.2.3. 

6.2.2  Other  Progressions 

The  feed  network  was  adjusted  for  0°  and  130° 
phase  progressions.  Figure  31  shows  the  result  of  0° 
progression.  Because  of  strong  coupling,  both  the  ar¬ 
ray  and  center  element  patterns  deviate  from  being  iso¬ 
tropic.  Figure  32  shows  the  result  of  180°  progression 
The  array  pattern  shows  four  nulls ,  as  depicted  by  the 
theory . 

6.2.3  Elevated  Center  Element 

In  order  to  observe  the  effect  of  higher  isolation 
the  center  element  was  elevated  from  the  array  plane  by 
a  quarter  wavelength.  The  isolation  of  the  90°  progres¬ 
sion  was  found  to  have  improved  by  3  db,  and  radiation 
patterns  were  plotted  in  Figure  33 .  An  improvement  is 
seen  in  the  center  element  pattern  as  expected  and  as 
explained  in  section  6.2.1. 

6.2.4  Effect  of  Other  Obstacles 

In  order  to  study  the  effect  of  obstacles,  loaded 
and  short-circuited  antennas ,  and  a  piece  of  straight 
wire  were  used  in  succession  to  simulate  various  situa¬ 
tions.  Figures  34  through  36  show  the  radiation  pat¬ 
terns  of  a  single  coaxial  dipole  in  the  neighborhood  of 
a  loaded  dipole,  a  piece  of  straight  wire,  and  a  short- 
circuited  dipole  respectively.  Pattern  distortion  is 
seen  to  become  worse.  As  a  comparison,  Figure  37  shows 
the  pattern  of  a  90°  progression  array  one  wavelength 
away  from  a  lr  .ded  dipole,  measured  from  the  array  edge 
Distortion  is  negligible.  The  reason  is  partly  because 
the  spacing  is  larger,  but  the  other  reason  could  be  be¬ 
cause  the  array  itself  is  0.6  wavelength  in  diameter, 
thus  the  effect  of  the  obstacle  tends  to  be  distributed 


among  array  elements ,  rather  than  concentrated  in  one 
element . 

6.3  Isolation 

For  the  S0°  phase  progression,  a  series  of  iso¬ 
lation  data  points  were  taken  and  mapped  onto  the  array 
plane.  Using  linear  extrapolation  between  points,  an 
equal-isolation  contour  plot  is  generated.  Using  the 
extrapolated  null  point  as  the  origin  and  the  lines  con¬ 
necting  it  to  antenna  locations  as  the  axis ,  displace¬ 
ment  away  from  the  null  was  marked  off  along  the  axis. 
Isolation  values  read  from  this  plot  compares  closely 
with  theoretical  results,  especially  the  trend  of  in¬ 
creasing  isolation  as  the  position  error  is  decreased. 
The  plot  is  shown  in  Figure  38. 

The  difference  between  theory  and  experiment  is 
largely  due  to: 

1.  The  theory  assumes  symmetrical  dipoles,  where¬ 
as  the  experiment  used  coaxial  dipoles. 

2.  Network  loss  is  present. 

3.  The  effect  of  an  absorber  ground  plane  which 
is  placed  one  wavelength  away  from  the  bottom 
end  of  the  dipole  to  minimize  the  effect  of 
supporting  structures  may  be  present, 

6,3.1  Value  of  High  Isolation 

A  high  isolation  of  55  db  was  measured  for  the  90° 
progression.  Higher  values  can  be  achieved  if  mechani¬ 
cal  tolerances  can  be  held. 

There  are  two  types  of  isolations  that  exist.  One 
is  the  isolation  between  the  array  and  the  center  ele¬ 
ment.  The  other  is  the  isolation  between  one  array  pro¬ 
gression  and  the  other.  The  first  is  controlled  prima¬ 
rily  by  position  precision.  Other  factors  that  do  enter 
into  the  consideration  are  array  phase  precision  and  ar¬ 
ray  amplitude  precision. 

Isolation  between  progressions  is  determined  mainly 
by  the  feed  network  matrix's  isolation.  At  present  the 
best  estimate  of  achievable  matrix  isolation  is  about 
30  db  for  a  transmission- line  network  and  40  db  for  a 
waveguide  network , 
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6.3.2  Isolation  Between  Two  Kalf-Wave  Dipoles 

In  existing  FAA  facilities ,  the  model  of  ‘two  di¬ 
poles  can  adequately  depict  the  isolation  between  an¬ 
tennas  .  Figure  39  shows  the  theoretical  plot  of  dipole 
isolation  as  a  function  of  spacing.  It  is  seen  that  in 
order  to  get  30  db  isolation,  the  dipoles  have  to  be 
spaced  4  wavelengths  apart,  and  to  get  55  db  isolation, 
the  two  dipoles  cannot  be  mounted  on  the  same  FM  tower 
or  similar  structures  of  that  size.  In  contrast  to 
that,  the  circular  array  needs  no  more  room  than  a  cir¬ 
cular  aperture  0.6  wavelengths  in  diameter. 

6 . 4  Gain 

The  circular  array  gain  at  90°  progression  was 
measured.  The  measurement  made  use  of  the  insertion 
loss  technique  and  Friis  transmission  formular.  The 
gain  was  found  to  be  0.5  db  above  an  omnidirectional 
antenna.  This  gain  value  includes  the  feed  network 
loss,  and  the  coupling  loss  to  the  center  element. 

6 . 5  Other  Capabilities 

Because  the  array  is  symmetrical  about  the  ver¬ 
tical  axis,  uniform  scanning  can  be  expected  if  fu¬ 
ture  needs  require  such  capability.  Beam  forming,  and 
simultaneous  beam  scanning  are  all  possible. 


7 1  CONCLUSIONS 


It  is  essential  to  improve  antenna  isolation  in  order 
to  up-grade  present  FAA  communication  systems.  The  concept 
of  using  progressively  phased  circular  arrays  to  achieve 
high  antenna  isolation  has  been  investigated  theoretically 
and  experimentally.  It  is  found  to  be  technically  feasible. 
Isolation  as  high  as  55  db  was  achieved.  Pattern  distortion 
is  not  severe,  as  compared  to  distortions  caused  by  nearby 
obstacles  in  existing  systems .  Space  needed  to  mount  the 
system  is  one  order  of  magnitude  less  than  what  is  required 
at  present  for  the  same  isolation. 

Disadvantages  as  compared  to  using  single  dipoles  are 
found  to  be  that  isolation  depends  on  position  accuracy. 

That  problem  can  be  solved  by  building  the  array  as  a  mod¬ 
ule.  If  weather  resistance  is  required,  a  radome  can  be 
used.  Because  of  the  feeding  network,  the  array  is  lossier. 
The  network  loss  is  estimated  to  be  in  the  order  of  1  db. 

Other  possible  disadvantages  are  that  the  array  could 
cost  more  than  single  antennas  because  of  the  additional 
components  in  a  phased  array;  and  an  increase  in  parts  may 
also  decrease  the  reliability  figure.  The  higher  manufac¬ 
turing  cost  is  there,  and  is  expected  because  the  array  can 
do  more.  Lower  reliability  figure  may  not  be  significant 
because  The  array  system  uses  only  linear  passive  components 

On  the  other  hand,  other  advantages  are  that  the  anten¬ 
na  system  is  complete  in  itself.  It  can  be  expanded  by  ei¬ 
ther  going  to  a  larger  matrix  network  [3]  or  using  the  modu¬ 
lar  construction  concept  and  stack  the  array  up  vertically. 
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RECOMMENDATIONS 


This  initial  study  has  proven  the  feasibility  of  the 
concept  of  progressively  phased  circular  array  for  antenna 
isolation.  It  is  one  of  the  few  that  are  applicable  to 
FAA  communications.. 

It  is  thus  strongly  recommended  that  the  study  be  ex¬ 
tended  to  investigate  the  limitations  of  matrix  feed  net¬ 
works,  and  effects  of  these  limitations.  Those  limitations 
that  should  be  investigated  are  loss.,  amplitude  balance, 
phase  precision,  and  isolation  between  phase  progression 
modes . 
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APPENDIX  A 


GENERALIZED  ADMITTANCE  MATRIX 

This  Fortran  IV  program  calculates  the  generalized 
admittance  matrix,  and  associated  admittance  and  imped¬ 
ances  . 


Plots  of  the  current  distributions  along  dipoles 
can  be  obtained  for  all  phase  sequences.  Necessary  sub¬ 
routines  are  included. 

INPUT  DATA: 

EH  -  half  height  of  dipoles  (height  of  monopole) . 
NN  -  number  of  dipoles . 

IPP  -  number  of  match  points . 

X(I)  -  X-coordinates . 

YU)  -  Y-coordinates . 

A  -  radius  of  dipole. 

I  -  index  of  antenna  number. 


OUTPUT  DATA: 

mum  M 

X(I),  Y(I) 

A 

CS(K.L) 

CSSIMUJ,KJ) 


defined  above, 
radius  of  dipole. 

the  generalized  admittance  matrix, 
short-circuited  admittances  and 
open-circuited  impedances , 


Assumptions  are  that  the  elements  are  thin,  center- 
fed  symmetrical  dipoles.  Surface  resistance  is  negligi¬ 
ble.  Kernel  is  approximated.  Moment  method  is  used. 


Some  of  the  subroutines  are  the  work  of  staff  and 
students  of  University  of  Mississippi. 
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IMPLICIT  COMPLEX (C ) 

COMPLEX  CSN(40,40) 

COMPLEX  OETERM, COEF ( 60 ) ,CV ( 64  J 
C0MPL6X*16  CS(40,40) 

COMPLEX  CSS JM( 10, 10} 

DIMENSION  AM AG  I ( 240 ) 

COMMON  X(5),Y{5).Z»A,AK,N»I 

C  A  IS  RADIUS,  EA  IS  HALF  HEIGHT  OF  CENTER  FED  DIPOLE, 

IPP  IS  NUMBER  OF  MATCH  POINTS*  NN  IS  NUMBER  OF  DIPOLES, 

M  +  N  ARE  FIELD  POINTS,  M=MATCH  PT,  N=D!POLE 
I  +  IP  ARE  SOIRCE  PTS,  I=DIPOL£,  IP=BASIS  SET. 

IF  THE  END  MATCH  POINT  SERVES  AS  A  SOURCE, IT  IS  ZERO,  BECAUSE 
OF  THE  BOUNDARY  COND l T ION ' WH ICH  STATES  THAT  THE  CURRENT  AT  THE 
END  OF  THE  a  IRE  IS  ZERO.  BUT  IF  THE  END  HATCH  PT  SERVES  AS  A 
FIELD  PCINT,  SINCE  IT  IS  NUST  A  LOCATION  FOR  COLLECTION  OF  MAG. 
VECTOR  POTENTIAL,  IT  IS  NOT  NECESSARILU  ZERO. 

EXTERNAL  CKS 
MT 1 =25 

AK=6.2H3L853 

C  READ  IN  EH  , NN,  IPP. 

READ(5,210)cH,NN, IPP 

210  FORMAT ( F5. 3 , 212 ) 

NNL1=NN-1 
ANNL  1=NNL 1 
ANNN=NN 
MM= I PP+ I 
KT=NN*MM 
MLNN=MT-NN 

C  REAC  IN  COORDINATES  CF  DIPOLES. 

READ { 5 , 7 ) ( X {  I), Y(I), 1  =  1, 5) 

7  FORMAT ( 2F8 .4 ) 

C  READ  IN  RADIUS  OF  ARRAY 

RE AO ( 5 , 2 ) A 

2  FORMAT ( F8.6 ) 

C  WRITE  OUT  CCOROINATES  OF  DIPOLE. 

rtR  I T  E  I  6 1  5 ) 

5  FORMAT! IOX, 'COORDINATES  OF  DIPOLE') 

LO  9  1=1, NN 

9  wR  I  TE ( 6 ,4 )  l  ♦  XU  )  » I  ,  Y  (I  ) 

4  FORMAT { 5X  ,  •  X ' , I  1 ,  •  =  « , F 10 . 4 , 5X , * Y ' ,  1 1  ,  •  =  » , F 10 . A ) 

WRITE (6,3) A 

3  F0RMAT15X, 'RADIUS  IS=',F9.6) 

ANN=  IPP 

UZ  =  EH/A.NN 
li=DZ/A*4. 

C 

C  INITIALIZE  CS(M,N) 


A- 2 


f.O  AC-0  .v  =  itM 
CO  AGO  N=l,MT 
A00  CS(V,NJ=CKPLX(0.,0. } 

CALCULATE  MATRIX  ELEM  BEGIN  WITH  THE  FIRST  HALF  SECTION. 


CO  101  N=1,NN 

IP  =  1 

ZN=0. 

ZNP l=0Z/2 • 

00  10A  M=i,MM 
AM=V 

Z- ( AM-1 . )*DZ 
DO  105  1=1, NN 
L  =  IP  +  (  I-D^IPP 
J=M+{N-i)*M.v 

CALL  CWEOF (CKS» ZN,ZNP1 , 1 1 ,CETERM) 
CS( J,L)=DETt«M 
105  CONTINUE 
10A  CONTINUE 


CALCU.ATE  FULL  MATCH  SECTION — MATRIX  ELEM.,  WHERE  I  AT  IS  MADE  0 


l’O  101  I  P=2  » IPP 
AP  =  I P 

ZNsQZ* ( AP- 1.5) 
ZNP 1=ZN+DZ 
1)0  L03  N  =  1 ,  M  M 
AM  =  V 

Z=( AM-1. )*UZ 


1)0  102  I  =  i  » NN 


L  =  I  P  +  (  I  - 1 )  <=  I  PP 
j=y  +  {n-  l ) 

CALL  CwEOF (CKS,ZN,ZNP1, 1 1 , DETERM) 
CS( J,L)=OETLRM 

102  CONTINUE 

103  CONTINUE 
101  CONTINUE 


CALCULATE  TERMS  CuE  TC  COSINES  (  EVEN  SYMMETRY) 


00  200  M=I,MM 
AK«M 

ZM= ( AM- 1 • )*CZ 
00  200  N=1,NN 
NINDX=M+(N-I  )*MM 
RS  =  COS(  AK*ZH/30. 

20C  CSl\lNDX»NNMPP+N)=CMPLX(0.  ,RS) 


, _  -*vv‘'^p 


o  rs  n  nnn 


I-MVrKT  MATRIX 


C 

CALL  CDMINl(M7,CS,40,CDeTM) 

V,.RIT£(6»989)CD£TM 

rtK17E(6,990) ( {CS{K,L),L=l,HT),K=l,MTl) 

L‘0  205  K=i.*MTl 
00  205  L=1»MT 
205  CSN(K,L)=CS(K»L) 

rt^ire(6,lMX(I),  1  =  1,5) 

W»ITC{6,1) C  y ( I J  *  1  =  1,5) 

1  FORMAT (5(2X,E12.6)  ) 

V.R  I  TEC  7, 1005)  (  (CSN(K,L)  »L=1,M7)  ,K=l,MTl) 

1005  FORMAT ( 20A4 ) 

CALCULATE  THE  DRIVING  COLUMN  MATRIX  FOR  THE  FIRST  ELEM 

DO  300  M=1,MM 
AM=M 

Z=  {  AM-1 .  )  <-'CZ 
RV=-S IN { AK*Z ) /60 . 

300  CV{M)=CMPLX(0. ,RV) 

CALCULATE  SHCRT  CIRCUIT  ADMITTANCES 

00  701  KJ  =  1,N.\ 

JK=K J*MM-MM 
DO  701  I J  =  1 » NN 
J=( I J-1)*IPP+1 
COEF( J)=CMPLX(0. ,0. ) 

CO  700  L=1,MM 
CHANGE =CS ( J  »  JK  +  L ) 

700  COFFI J)=CGEF( J)+CHANGE  *CV(L) 

701  CSSIMI  IJ,KJ)=CCEF{ J) 

WRITE (6, 921) ( (I J.KJ.CSSIMI I J , K J ) , K J=l ,NN ) » I J= l» NN ) 

C 

C  C.ALCULAlt  OPEN  CIUCUIT  IMPEDANCES 
C 

CALL  CMINI(NN,CSSIM, 10, DETERM) 

„RITE  (6,991  JDFrERM 

WRITE  (6, 920  I  (IJ,KJ,  CSSIMI  I J  ,KJ  )  ,KJ  =  l ,  NN ) ,  IJ=  1 ,  NN ) 

920  FORMAT ( 4 ( '  Z  •  , 2  1  1  ,  «  =»,2E12.5)) 

921  FORMA  r ( A ( *  Y*,2I1,'  =  ',2E 12.5)) 

C 

C  °EbIN  PHASE  SEQUENCES. 

C  JJ=i  FOR  0,2  FOR  90,3  FGR  180, A  FOR  270. 

00  500  JJ= 1 ,NNL 1 

r* 

C  DRIVING  MATRIX  FOR  SUBSEQUENT  ELEMENTS. 

C. 


A~4 


JJ J=A JJ 
t'.O  301  M=1,MM 
C 

C  OR  I V IMG  6LEMS  OF  CENTER  ANTENNA  IS  THE  SAME  AS  ANT  =1. 

C 

CV(NNl1*MK+M)=CV(M) 
on  301  N=2  t.N'MLl 
NLl=N-l 
ANU=NlL 

A.MGLE=AK/ANNL1*AJJ«ANL1 
301  CV  (NL1*MM-*M  )  =CV  ( H )*CEXP  (C^PLXtO.  »ANGLE ) ) 

■NCENT-NNL1*  1 PP+ 1 
CGcFll)=CMPLX{0.,0.) 

CGEF { NCENT } =CCEF ( 1 ) 

MTLMM=MT-MM 

C 

C  FEtC  POINT  I  AT  ELEK.  =1  +  THE  CENTER  ELEH.  HAVE  NO  DRIVE  FROM  CENT . 
r 

Of)  331  L  =  1  » MTLMM 
CHANGE=CS( l,L) 

COEF { 1 ) =COEF l 1 J+CHANGE  *CV(L> 

CHANG2=CS ( NCENT ♦ L ) 

331  CO EF(NCENT)=CCEF( NCENT) +CHANG2  *CV(L) 

AMPL l =CAriS ( COEF ( 1 )  ) 
wRITE{6,9lO)CGEF(l) ,AMPLI 

AMPL I -CABS ( COEF ( NCENT  )  } 

WRITE (6,911) COEF (NCENT) , AMPL  I 
CZ  L2=C(JEF { NCENT) 

MTUMTLMMU 

C0FF(NCFNT)*CMPLX(0.,0. ) 

c 

C  F 6 E n  PT.  I  AT  THE  CENTER  ELEMENT  WITH  NO  DRIVE  FROM  CIRC.  ARRAY. 

C 

DC  3  32  L  =  M  T 1 , M  T 

332  COEF  (NClNTXCEF  (NCENT )+CS(NCENT,L)*CV(L) 

YL~ . 0 1 

C V?  =  -C7. 1 2/ (  CHEF  (NCENT  )+YL) 

Ampl l -CABS (COEF (NCENT)  ) 
kRITEI 0,912 )CCEF (NCENT) , AMPL I 
AMPL  I =C ABS ( C V2  ) 
wRITE(6,9l3)Cv2fAMPLI»YL 
JO  333  L=MT1,MT 

333  CVlL )=CV(L)*CV2 
C 

C  CALCULATE  CCFFFICIENTS 
r 

DO  3  20  J  =  I »  VLNN 
COEF( J)=CMPLX(0. »0.  ) 


A~5 


o  o  o 


DO  321  L=1  »MT 
CHANGF=CS! J,L> 

321  COEF  ( J  )  =CGEF  ( J  )  +CFiA\‘GE  *CV!L) 

AM AG  I !J)=CABS!COFF(J)  ) 

AM AG  I (  MLNN+J ) =  REAL ( COEF { J) } 

AMAG  I  { 2*NLNN+ J  )  =  A I  PAG  ( COEF  (  J  )  ) 

4f'.AGU3*l*LNN+Jl=0. 

32C  WRITE!6,999)CC£F! J)  ,AMAGI<  J) 

CALC  INPUT  IMPEDANCES  AND  ADMITTANCES 

DO  600  1=1, NN 
AlLl=l-l 

ANGLb  =  AK/AN.\Ll*AJJ*AILl 
MMI=MM#I-MM+1 
1 1  PL  I  =  ( I-i)*lPP 

LETtRM=COEF { 1 1  PL  1  +  1 ) /CEXP ICMPLX ( 0» , ANGLE ) ) 

4MPLI=CABS!0ETERM) 

WRITE (6, 904)1 ,CFTERM,AMPLI 
DETERM=CMP(.X{  1.  ,0.) /DETERM 
AMPlI=1./AMPLI 

600  WRITEI6, 903)1, DETERM, AMPLI 

SOO  CALL  RPLOT ( AMAGI ,4,HLNN,i» , 1 • , J  J J ) 

900  FORMAT  l  *  1=  ',13,*  CS(*,2I3,»)  =  S2E11.3) 

901  FORMAT { *  IM  =  *,13,'  Z=  *,£11.3) 

902  FORMAT { 1  I P , LK , ZNP 1 , 1 1 , N , =  * , 1 3 , 2E 1 1 . 3 , 2 14 ) 

903  FORMAT  (  *  ZP,I2,')  =',3E14.8) 

9C4  FORMAT  (  '  Y{',12,’)  =',3E14.8) 

910  FORMAT! •  Yll  =‘,2E14.8,‘  MAGNITUDE  =»»E14.8) 

911  FORMAT ( '  Y12  =',2E14.8,»  MAGNITUDE  =',E14.Q) 

912  FORMAT ( '  Y22  =',2E14.8,»  MAGNITUDE  ='?E14.8) 

913  FORMAT!'  V  =',3614.8,*  Y  LOAD  =*»E14.8) 

989  FORM AT ( 4fcH  INVERTED  MATRIX  WHCSE  NORMALIZED  DETERMINANT  ,2016.6) 
99C  FORMAT! 10D11. 3) 

991  FORMAT { 46H  INVERTED  MATRIX  «HCSE  NORMALIZED  DETERMINANT  , 2E16. 6 ) 
999  FORMAT { 3E16.6) 

1000  CONTINUE 
CALL  EXIT 
END 

COMPLEX  FUKCT I CN  CKS(ZP) 

COMMON  X  (  5  )  ,  Y  { !j ) , Z , A , AK,N, I 
PR=A**2+(X( I )-X(N) )**2+!Y! I )-Y!N) )**2 
Rl=S0RT{RR+(Z+ZP)**2) 

R?=SQRT (RR+(Z-ZP)*#2) 

RK=COS( AK*R15/Rl+CCS(AK*R2)/R2 
AIK--SIN! AK*Rl )/Rl-SIN!AK*R2)/R2 
CKS=CMPLX(RK,AIK) 

RETURN 

END 


A- 6 


.UBRCUTINE  CrtEOF  (CF»XL,XU,NX,CANS)  CWEOFOOl 

* 

CWEDF  IS  A  SUBROUTINE  WHICH  WILL  NUMERICALLY  INTEGRATE  A  USER  * 

SUPPLIED  FUNCTION  BETWEEN  SPECIFIED  LIMITS.  (SINGLE  PRECISION)  * 

CF  -  NAME  OF  FUNCTIC .  SUBPROGRAM.  MUST  BE  LISTED  IN  AN  * 

EXTERNAL  STATEMENT.  * 

XL  -  LOWER  LIMIT  OF  INTEGRATION  * 

XU  -  UPPER  LIMIT  OF  INTEGRATION  * 

NX  -  APPROXIMATE  NUMBER  CF  NODES  AT  WHICH  TO  EVALUATE  FUNCTION  * 


CAvS  -  RESULT  OF  INTEGRATION 

PREWARED  BY  MICHAFL  G.  HARRISON 

IMPLICIT  CC-MPLEX*a  (C5 
RFAL*d  CXDX t XX 

RCAL*4  CW16J/82. ,216. ,27. ,272., 27. ,216. 

IF(NX.Lt.G)  GC  TO  900 

N=< (\X+4)/6)*fc+l 

IiX  —  (  XU-XL  )  /FLOAT  { M—  1 ) 

CXCX=OdLE(CX) 

•Nw  I  X=N/6 
X  =  XL 

CA.\IS=-CF(X)*41.0 
on  see  mx= i , nw i x 

1  0  700  r » X  =  1  f 6 
CANS=CAi\S  +  Cw(KX)*CF(X) 

XX  =  OBL  E ( X ) 

X=5NGL ( XX+CXDX ) 

CONT I NUc 
CONTINUE 

CANSs(CANS  +  4l.0*CF(X)  )$0X/ 140.0 
RETURN 

wP I T  E i 6  *  90 1 )  N 

rURMAT ( 1 OERROR  IN  CALLING  PARAMETER  **** 

RETURN 

ENC 

SUBROUTINE  RWLCT 


E.E.  DEPT  JUNE  22,  1972  * 

tx 

CWE0F0G2 
CWEDF0Q3  ; 

CWEDF004 
CWEDF005 
CWEDF006 
CWEDF007 
CWEDF008 
CWEDF009 
CwEDFOlC 
CWEDF01 1 
CWEDF012 
CWEDF013 
CWEDF014 
CW6DF015 
CWEDF016 
CWEDF017 
CWEDF018 
CWEDF019 
CWEDF020  } 
CWEDF021  { 

N  =  • » 1 5  t '  *****«//)  CWEDF022 

CWEDF023 
CWEDF024 


PURPOSE 

PLOT  SEVERAL  DEPENDENT  VARIABLES  AGAINST  ONE  INDEPENDENT 
VARIABLE 


USAGE 

CALL  RPLCT  (A,NC»NR»X1,XI,IN0) 

INSCRIPTION  OF  PARAMETERS 

A  -  THE  ARRAY  TC  BE  PLOTTED.  EACH  COLUMN  CONTAINS  A 
VARIABLE  TO  BE  PLOTTED 

-  THE  NUMBER  CF  COLUMNS  IN  A 

-  THE  NUMBER  CF  ROWS  IN  A 


NC 

NR 


XI  -  THE  FIRST  VALUE  OF  THE  INDEPENDENT  VARIABLE 
XI  -  INCREMENT  CF  THE  INDEPENDENT  VARIA8LE 
INO  -  CHART  DUMBER  (3BIGITS  MAXIMUM) 


REMARKS 

NONE 

SUPROUflNES  AND  FUNCTION  SUBPROGRAMS  REQUIPf.C 
NONE 


SUBROUTINE  RPLCT  (  A ,NC ,NP. , XXI , XI ,  INO) 

REAL  MlN,MAX,A(l) 

INTEGER  IR( 12) »LINE( 101 ) ,  BLX 
Xl-XXI 

K)RVAT  ( 18H 1  CHART  NUMBER  ,  1 3 » / ,  IHO,  E  Lr 6 *7 IX  ,£  15 . 6 , 13X ,  7HX  VALU 

c,//) 

FORMAT  {  i  H  ,101Ai,8X,EL5.6) 

FORMAT!  IH0,80X, 16, 15F  POINTS  PLOTTED) 

INITIALIZE  VARIABLES 


ELK  =2**  30 

IR(l)=2**30+2**27+2**25+2**24  PLOT 

IR(2)=2-:*30+2**29+2**27+2**25+2**24  PLOT 

I P,  (  3)=2**3C+2**27+2**2fc+2**25  PLOT 

IR (41=2**30+2**28+2**27+2**26  PLOT  •*' 

IR(5)=2** 30+2**27+2**25+2** 24  PLOT  *  L  * 

IR ( 6 ) =2** 30+2**29+2* *27+2**25+2**24  PLOT  '2* 

IR(7)=2**30+2**27+2**26+2**25  PLOT  *3 

H (91=2**30+2**28+2**27+2**26  PLOT  ’4' 

IK (91=2**30+2**27+2**25+2**24  PLOT  *5' 

1R(  101=2**30  +  2**29  +  2**27+2**25  +  2**24  PLOT  ‘6’ 

IK(  ll)=2*v3C+2**27*2**26+2**25  PLOT  '7' 

IR( 12)=2**30+2**?b+2**27+2**26  PLOT  '8* 


00  4  1=1,101 
LINE! I  )=BLX 

continue 

LOCATE  MIN  ANC  M AX  VALUES 


N=NC*NR 
M l N= A { 1 ) 
MAX  =  A (  l ) 


A- 8 


on  e  1=2  »n 

I F ( A { I ) — K IN)  5,6,6 

5  MIN  =  A { I ) 

6  CONTINUE 

IF ( A ( I )— MAX )  8,0,7 

7  KAX=A(I) 

0  CONTINUE 

C 

C  FOR  SINGLE  VALLEC  ARRAY  SET  LIMITS 

C 

IF l  AX-MiN)  9,9,  IG 

9  MAX=A(l)+i.O 
KIN=A ( 1 )  —  1 «  0 

10  CONTINUE 
C 

WR I TE (6 , 1 )  INC, MIN, MAX 
C 

C  BEGIN  PLOT  LOCP 

C 

00  15  1=1, NR 
DO  14  K=1,NC 
1FIK-1)  11,11,12 

11  KSA=I 

GO  TG  13 

12  KSA=KSA+NR 

13  CONTINUE 

KPNT=( A (KSA)-M IN )/{ MAX-MI N) #100.0+1.5 
LINE!  KPNT )  =  IR ( K ) 

14  CONTINUE 

wft I T  £ ( 6 , 2 )  LINE, XI 
X1=X1+XI 
00  15  L  =  1 , 1 0 1 
LINE l L ) =BLK 

15  CONTINUC 
r 

rtR  I  TF { 6 , 3 )  NR 
C 

kE  TURN 
FNC 

SUBROUTINE  CMINl  ( N , A ,\D I M, GETERM )  CMN10010 

C  * 

C  C MINI  IS  A  SUBkCUTINE  WHICH  WILL  ACCEPT  A  SINGLE  PRECISION  COMPLEX  * 

C  MATRIX  AND  R  E I  URN  THE  INVERSE  OF  THE  MATRIX  IN  ITS  PLACE.  THE  * 

C  SUBROUTINE  WILL  ALSO  COMPUTE  THE  NORMALIZED  DETERMINANT  OF  THE  MATRIX.* 

v.  't  *  THE  CROER  OF  THE  MATRIX  TO  BE  INVERTED  * 

C  A  COMPLEX  DOUBLE  PRECISION  INPUT  MATRIX  <  DESTROYED )  * 

C  THE  INVERSE  OF  A  IS  RETURNED  IN  ITS  PLACE.  * 

C  ND I M  -  THE  SIZE  TO  WHICH  A  IS  DIMENSIONED  IN  THE  CALLING  PROGRAM  * 

C  OEfcrtM  .'mE  NORMALIZED  DETERMINANT  WHICH  IS  CALCULATED  BY  THE  * 


A-9 


o  o  o  n  o  n  o  n  n  o  o  o 


DETERM  -  THE  NORMAL  I  ZED  DETERMINANT  OF  A  WHICH  IS  RETURNED  * 

PREPARED  BY  MICHAEL  G.  HARRISON  E«E*  DEPT  JUNE  23,  1972  * 

* 


COMPLEX  A!NDlN»NOIM) , PIVOT! 100) ,AMAX,T, SWAP, DETERM, U 
INTEGERS  i PIvCT ( 100  ) ,  INDEX  U00 , 2 *- 
REAL  TEMP, ALPHA! 100) 

INITIALIZATION 

OETERM.  =  CMPLXtl.  0,0.0) 

00  20  J=l,\ 

ALPHA! J)=O.CDC 
00  10  1=1, N 

10  ALPHA! J)=ALPHA!J)+A!J, I )*  CON JG ( A ( J , I ) ) 

ALPHA ( J )  =  SGRT ! ALPHA ( J ) ) 

20  I PIVCT { J ) =0 
CO  600  1=1 , N 

SEARCH  FOR  PIVCT  ELEMENT 

AMAX=CMPLX!C. 0,0.0) 

DO  LOb  J  =  1 .  K 

IF  (IPiVOT(J)-l)  60,  1C5,  60 
60  00  100  K  =  1 , N 

IF  (IPIVOT(W)-I)  80,  LOO,  740 
80  TEMP=AMAX*  CON JG < mMAX ) -A { J , K ) *  CGNJG! A! J,K) ) 

I F ( TEMP ) 85 , 85, 100 
85  IRGw=J 
ICOLUM=K 
AMAX  =  A! J,K  ) 

100  CONTINUE 
105  CONTINUE 

IPIVOT! ICCLUM)=IPIVCT( IC0LUM)+1 

INTERCHANGE  RCHS  TO  PUT  PIVOT  ELEMENT  CN  DIAGONAL 

IF  {  iROVi-ICCLUM)  1  AO ,  260,  140 
140  DETERM=-DETERM 
DO  200  L=l,N 
SWAP=A ( 1  ROW  ,  L  ) 

A  (  I  ROW ,  L  )  =  A  (  ICCLl'M,  L  ) 

200  AUCOLUM.L  )=SV<AP 
SWAP=ALPHA( IRGW) 

ALPHA!  IxOW)=At.PHAl  ICCLUM) 

ALPHA l I CCLUM. ) =SWAP 
2fcO  INDEX! I , 1)=IRCK 
INCEX ( I ,2) = ICCLUM 
PIVOT!  I  )=A(  1CC-LUM, ICCLUM) 

U  =  PIVCT! I  ) 


CMN10020 

CHN10030 

CHN10040 

CMN10050 

CHN10060 

CMN10070 

CMN10C8C 

CMN10C90 

CNN10100 

CMN10110 

CMN10120  * 

CHN10130 

CMN10140 

CMN10150 

CMN10160 

CMN10170 

CMN101P.0 

CMN1G1 90 

CMN102G0 

CHN10210 

CMN10220 

CMN10230 

CHN1G240 

CMN10250 

CMN10260 

CMN10270 

CMN10280 

CMN10290 

CMN10300 

CMNL03I0 

CMN10320 

CMN10330 

CHN10340 

CMN10350 

CMN10360 

CMN10370 

CMN10380 

CMN10390 

CMN10400 

CMN1G410 

CHN10420 

CMN10430  . 

CMN10440 

CMN10450 

CMN10460 

CMN10470  • 


A-10 


r  •  r>  o  o  r>  o  <>  o  no  on 


. STcRM  =  nETE<vM*L 

O.-  TERM =GET£RM/ AlPHA  ( ICCLUK ) 

T£PP=PI vOT 1 1  )*  CONJGIPIVCTI I } } 

Ir(rs:KPJB3C,720,330 

MViuE  PIVCT  RC.W  f)Y  PIVOT  ELEMENT 


33C  - ( ICGLUM » ICCLUM )  =  CKPLXI1. 0,0*0) 
»<•*  35C  L=i  t »\ 

o  =  prvcrn  i 

350  At  ICCLuM  « L )  =  A  (  ICCLUM, L I  /U 
rfcCLCE  NON-PIVOT  KUWS 


80  HO  550  Ll=i,N 

IFIL1-ICQLLM)  400,  550,  40C 
40C  T=A{L1, ICCLUM) 

Art  LI ,  1C0LUH )=  CMPIXIO. 0,0.0) 
DC  450  L  =  I » \ 
t.  =  A{  ICULLPjL) 

450  A { Ll , L )  =  A ( L I , L ) — U*T 
55C  COMINUS 
0  CONTINUfc 


INTERCHANGE  COLUMNS 


20  DO  710  1=1,. 

L=N+L-T 

IF  { INDEX ( L , 1 )- INDEX ( L,2 ) )  630,  710,  630 
630  JRCW=INCEX(1,1 ) 

JC0LUK= I NOFX I L , 2 ) 

00  705  K  =  1,N 
SwAP=A(K,jROX) 

■*  (K»JRuW)  =  A(K,JCOLUv) 

A  ( K , JCGLUM ) =SW AP 
705  CONTINUE 
71C  CONTINUE 
RETURN 

720  WRITE  (6, 730 

73C  FORMAT ( 20H  MATRIX  IS  SINGULAR) 

740  Rf  TURN 
E\C 

SuuROuTINE  CUM  INI  IN,A,N01M,0ETERM> 


CMN1C450 

CXMC-9G 

CKNICsCO 

CHNiOSK; 

CHN1C520 

CM.MC53G 

CRMO^O 

CHN1055C 

CHNlG“i. 

OMNIOo/ 

C^NlOSo." 

CMN106, T 

CMN106G:i 

CMNICoiC 

CMNlOtE  0 

CMN1G63G 

CHN1C6A0 

CHN10A*  ‘ 

CHNlCAcu 

CMN10670 

CHN1066G 

CHNi0fc5u 

CMN1G7CC 

CHN1371C 

CMN1072C 

CMN10730 

CHN1C74G 

CHN1C75U 

CMN10760 

CMN1077G 

CMN107e0 

CMN10790 

CHN 10800 

CMNlOciO 

CHN10M2C 

CHN1G83C 

CMN10e-'-C 

CHMGS50 

CMN10B60 

CMN10870 

CMN1086C 

CKN1039G 

COM  10010 


:CM!M  IS  A  SUBROUTINE  WHICH  WILL  ACCEPT  A  0C-U9LE  PRECISION  COHPLEX  * 
"ATR I X  AND  RETURN  THE  INVERSE  OF  THE  MATRIX  IN  ITS  PLACE.  THE  * 

SUKRuuT 1 NC  h I LL  ALSO  COMPUTE  THE  NORMALIZED  DETERMINANT  OF  THE  MATRIX." 
C  N  -  THE  CROER  CF  THE  MATRIX  TO  BE  INVERTED  * 

C  A  COMPLEX  DOUBLE  PRECISION  INPUT  MATRIX  (DESTROYED)  * 


A- 13. 


non  (u  non  r«  o  r>  o  r.  n  r- 


THfc  INVERSE *GF  A  IS  RETURNED  IN  ITS  PLACE. 
NCI"  -  THE  SIZE  TO  WHICH  A  IS  DIMENSIONED  IN  THE 
CfcitRP  -  THE  NORMALIZED  DETERMINANT  OF  A  WHICH  IS 
PREPARED  8Y  MICHAEL  G.  HARRISON  E.E.  DEPT  JUNE  23 

CCPPi.EJC*16  AINCIMfNDIMIjPlVCTUOO)  tAHAX,T,SWAPsDETERHtU 
C0MPLEX*16  CCKPLX,DCQNJG»CDINV,COXXXX 
1NTEG£R*4  I P IVCT ! ICO ) , INDEX (100,2) 

REALMS  ALPHA! 100), TEMP 


INITIALIZATION 

uETERK  =  OCMPLXI i.C*0 ,0.0+0) 

DO  20  J  =  l*N 
ALPHA! J ) =C.GOO 
DO  10  1=1, N 

10  ALPHA! J)=ALPHA(J)+  «  *» I ) *OCONJG! A i J , I ) ) 

ALPHA! J)=DSCRT( ALPHA! J) ) 

20  IPI VGT ! J ) =0 
DO  600  1=1, N 

SEARCH  FOR  PIVCT  ELEMENT 

AMA<  =  0CKPLX(C.0+0,C.0+0) 

CO  105  J=1.,N 

IF  IIP.IVOT!  J)-l)  60,  105,  60 
60  DO  100  K=l,\ 

IF  (IPIVCT(K)-I)  80,  100,  740 
t  TFMP=AMAX*DCONJG!AMAX)-A< J,K)*DCGNJG{A( J,K) ) 

IF  I  TEMP )85, 85,100 
83  IRCV>= J 
ICCLUM=K 
AMAX=A! J,K) 

IOC  CONTINUE 
105  CONTINUE 

IP  I VGT! ICCLLMJ  =  IPIVCT(ICCLUM)  +  1 

INTERCHANGE  RCRS  TC  PUT  PIVOT  ELEMENT  CN  DIAGONAL 

IF  (  IRGw-ICCLUM)  140,  260,  14C 
14C  DETERM=-D6  TERM 
CO  200  L=l,N 
SWAP= A ( IRCW.L) 

A!  IROw*U=A(  !  CCLUM ,  L  ) 

2CC  A(  ICOLU'v,L)=SwaP 
SwAP= ALPHA ( I  ROW) 

ALPHA!  IRCW)=ALPHA(  ICCLUM.) 

AlPHA(  ICGLUM)  --swap 
260  INDEX! !,l)=IRCfc 


CALLING  PROGRAH  * 
RETURNED  * 

,  1972  * 

$ 

C0K10Q20 
C0K1C030 
CDM10040 
CDK1G050 
CDM10060 
C0H1C07O 
CDM10080 
C0H10090 
COHIOIOO 
CDMI0110 
CDH10120 
CDM10130 
CDK10140 
CDM10150 
CDH10160 
C0M10170 
CDM10180 
CDK10150 
CDM10200 
CDH10210 
CDM10220 
CDH10230 
C0MI0240 
CDM10250 
CDM10260 
CDH10270 
CDM10280 
CDM1029C 
CDKI0300 
CDH10310 
CDMI0320 
CDM10330 
C0MI034O 
CDM10350 
CDH10360 
CDM10370 
CDM10380 
COM  10390 
CDM104C0 
CDM10410 
CDM10420 
CDMI0430 
CDM10440 
CDM10450 


r.  n  o »  nno  c»  r  >  r> 


ivc-.xc  t  ,?j-ir.rt  um 
piveti  d=a»  ici.lum,  icclum) 

U  =  FltfOT(i) 

DFTERM  =  GETERM4U 

C-E  Ii:«M=CETERH/ALPKA  UCOLUM ) 

TEMP=PlvflT  ( !  )*CCO.NJG(  PIVOT { I ) ) 

IF {TEMP) 330,720, 330 

DIVIDE  PIVCT  RCW  t>Y  PIVOT  ELEMENT 

330  A(  ICCLUM » ICCLUH)  =  CCMPLXll.D+O, 0.0+0) 

uO  L=  1 1  *\ 

o  --  :>ivt.rm 

350  A ( ICGLltM , L )  =  At I CGLUH, L ) *C0XXXX (U > 

SLcCuCE  NON-PI VCT  RCHS 

380  CO  65C  L 1  =  I , N 

I P ( L 1- ICOLLM )  400,  550,  400 
400  T=A(LI, ICOLLM) 

A  (  L  1 ,  ICGLUM )  =CCMPLX  (0.D+0,0.C+0) 

!’0  450  L  =  1,N 
li  =  A  UCOLUM, L) 

45G  A  { L I ,  L  )  =  A(L1,L)-L*T 
550  CONT I NUh 
0  CGNTINUt 

INTERCHANGE  COLUMNS 

620  CO  710  1  =  1, N 
L=N+  1-  I 

IF  { INDEX { L , 1 ) - I N06X ( L , 2 ) )  630,  710,  630 
630  JRrw=INOEX  f.,1) 

JCCLLMUNCEXU  ,2  ) 

00  706  K=  1  i N 
SWAP=A{K,  JRC'W) 

A { K , JRCW ) =A ( K , JCClUM ) 

A(K»JCOLU.m)=S’aAP 
705  CONTINUE 
71C  continue 
RETURN 

12  0  WRITE  (6, 730) 

73C  FORMAT ( 20H  MATRIX  IS  SINGULAR) 

740  RETURN 
END 

COMPLEX  FUNCTION  Ci)XXXX*16<  A) 

COMPLEX’8'16  A ,  CCMPLX 

REAl*H  AR,  A\,ARINV,AIINV,0A3S 

AR  =  A 


CDK10**60 
CDM1C470 
CDMI04LO 
C0H1049C 
COM  1 C500 
CCMICsLO 
CDM1G520 
COM 1C 530 
CDH1G540 
CDM10550 
CDH1C560 
COMirsvc 
C0MlC5cC 
CDM10590 
COM 10600 
COM 1 0610 
COM l 0620 
CDHlG6*.i> 
COMIC-, 40 
COMlOU-U 
CDMICP60 
C0M1O6/O 
c  oM  1  ’  * 

ca^‘.C6  , 
COmi;  a.< 
COM 10  / 1 0 
COHIO  I'iik 
COM  1 C 7  Hi 
C0H1074U 
CDM10750 
COM 10760 
C0M1077O 
CUt*  l  vhW 
CDMluH). 
f  t)Mlf ,  t  i 
i  ,.**  i  if  A  . 

(  i>#  i<>r< 
CUMh'H  h> 
cnnuiMi, 

l  '  4 

criMit  lf, 
Ll‘M  Hifi  .*•* 
COM  lie  .  i. 
COM  (1(1*14, 

Cl  *  ft*. o( 
C.  M 

C  uM »  > : ‘ 

H< 

COM  1 U 
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AI=-{0.0+Cfl.D+0)*A  C0K10950 
IF IQABSI AR )  . LC.  1 .0-30) AR=0.C+0  CDH10960 
IFIDABSiAp),  .LE.  i.C-30) AI=0.C+0  CDK10970 
ARINV=AR/(AR*AR+AI*AI)  COM 10980 
AIINV=-AI/(AR*AR+AI*AIJ  CDH10990 
COXXXX=OGKPLX( ARINVf AI IMV)  C0H11000 
RETLR.'J  C0HII01G 
EN'L  CDM11020 
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APPENDIX  B 


MOMENT  METHOD  DATA  REDUCTION  BY 
SEQUENCE  FUNCTIONS 


ABSTRACT 


The  generalized  admittance  matrix  of  wire  antennas 
obtained  by  applying  the  moment  method  is  usually  large, 
and  cumbersome  to  handle.  In  order  to  facilitate 
repeated  calculations  to  find  field  and  circuit  quantities 
from  the  matrix,  sequence  functions  are  used.  The  approach 
avoids  the  necessity  to  store  the  entxre  admittance  matrix 
and  also  reduces  subsequent  computing  efforts.  It 
supplies  current  distributions  and  field  patterns  of 
each  sequence  as  an  intermediate  step,  and  thus  provides 
otherwise  unobtainable  insight  into  the  performance  of  an 
antenna  system,  especially  that  of  a  conformal  array. 
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MOMENT  METHOD  DATA  REDUCTION 
BY  SEQUENCE  FUNCTIONS* 

Bing  Chiang,  Howard  University,  Washington,  D.C. 


In  applying  the  moment  method  for  solutions  of  wire 
antennas,  the  labor  is  in  calculating  the  impedance 
matrix,  and  in  matrix  inversion  to  obtain  the  admittance 
matrix.  Matrices  encountered  in  this  solution  method 
are  generally  very  large.  To  calculate  field  and 
circuit  quantities  directly  from  these  matrices  is  an 
equally  laborious  task  due  to  sheer  bulk  of  data.  In 
order  to  reduce  data  handling,  the  approach  of  using 
generalized  sequence  functions  is  suggested. 

Sequence  functions  have  been  used  by  King,  Mack  and 
Sandler £ to  analyse  symmetrical  circular  arrays. 
However,  antennas  and  arrays  often  have  arbitrary 
shapes  and  sizes  and  are  often  unsymmetrical ; 
therefore,  a  generalized  sequence  function  theory  is 
needed.  The  theory,  developed  below,  is  valid 
wherever  super-position  hulds. 

Current  distributions  on  wire  antennas  are  shown  in 
the  moment  method  as*-2-* 

[In]  »  [Ynm]  [vm]  CD 

where,  I  is  the  current  distribution  along  the  antenna, 

V  is  the  generalized  admittance  matrix, 

V  is  the  generalized  voltage, 

n  is  the  index  of  the  basis  set  which  goes  from 
1  to  N, 
and  lastly, 

m  is  the  index  of  the  testing  set  which  goes 
from  1  to  M. 

The  size  of  the  generalized  admittance  matrix  is  thus 
NxM,  where  N  j<  M. 

If  there  are  P  ports  in  the  antenna  system,  then  the 
matrix  size  can  be  shown  to  reduce  to  NxP  when  sequence 
function  is  used.  The  saving  is  thus  M/P  fold. 

Formulation  of  the  sequence  function  is  based  on  that 
any  applied  voltage  vp  in  a  set  of  P  voltages  can  be 
specified  by  a  linear  combination  of  P  sequence 
voltages : 

V  P-1  f  y 

P  =  2  Aw  exP  [j  2n(p-I)m/P]  (2) 

m  =  0 

where  is  the  complex  coefficient  of  the  mth 
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sequence.  Shelton- Butler  matrix  [33  is  the  parallel 
physical  system  that  performs  sequence  generation  and 
summation  as  shown  in  equation  (2).  An  array  fed  by 
such  a  matrix  is  shown  in  Figure  1.  In  this  figure, 
only  the  n/2  sequence  is  excited,  so  that  the 
progressive  phase  going  from  port  to  port  can  be 
illustrated  clearly.  The  random  orientation  of  antennas 
and  location  of  antenna  ports  depict  the  general 
applicability  of  the  theory.  When  all  matrix  input 
ports  are  excited  simultaneously,  the  voltage  at  the 
antenna  port  is  thus  the  Vp  shown  in  equation  (2).  Any 
set  of  Vp  can  be  completely  specified  by  a  set  of  A(m). 

Using  the  moment  method,  current  distributions  on 
antennas  can  be  calculated  from  the  generalized 
admittance  matrix.  Once  the  current  distribution  for 
each  voltage  sequence  is  calculated,  the  response  of 
any  arbitrary  set  of  voltages  vD  can  be  found  by 
scaling  and  superpositioning  .  And  because  vp  can  be 
synthesized  by  scaling  and  superpositioning,  it  is 
obvious  then,  the  only  information  worthy  of  storage  Is 
the  current  distributions  caused  by  a  set  of  normalized 
sequence  voltages.  There  are  N  data  points  for  each  set 
of  distributions,  and  there  are  P  sequences.  Total  data 
points  are  thus  equal  to  N  x  P,  which  is  a  reduction 
from  N  x  M. 

In  applying  sequence  voltages,  data  stored  are  not 
necessarily  limited  to  current  distributions.  If  one 
so  wishes,  he  may  store  the  complex  field  pattern  for 
each  mode  instead,  and  obtain  the  desired  field  pattern 
by  again  scaling  and  superpositioning. 

As  an  example,  a  circular  array  of  dipoles  as  shown  in 
Figure  2  with  an  array  radius  of  0.3  wavelength,  and  a 
passive  element  near  the  center  was  analyzed.  For 
simplicity  of  analysis  all  elements  were  made  equal. 

They  have  a  wire  diameter  of  0.025  wavelength.  The 
passive  element  was  displaced  from  the  array  center 
by  .03  wavelength,  and  was  loaded  by  a  100  ohm 
resistor.  Each  half  of  a  dipole  was  divided  into  5 
pulses,  with  end-sections  having  half  the  width. 

Field  and  phase  patterns  of  each  sequence  were  calcu¬ 
lated.  They  are  plotted  in  Figure  3.  Let  it  be 
assumed  that  the  desired  feed  voltages  are  those  shown 
in  Table  1.  Since  the  sequences  form  an  orthogonal  set, 
the  desired  sequence  voltages  can  be  found,  and  are 
listed  in  Table  2.  The  resultant  pattern  is  obtained 
by  scaling  and  superpositioning  the  patterns  shown  in 
figure  3  and  is  plotted  in  Figure  A. 
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APPENDIX  C 


MISCELANEOUS  COMPUTER  PROGRAMS 


INPUT  DATA: 


IPP 

NM 

EH 

CS(I,J,K) 

I 

J 

K 

X(L) 

Y(L) 

L 

RAD 

PERRX 

PERRY 


number  of  match  points . 
number  of  dipoles, 
half  height  of  dipoles, 
inverted  short-circuited  matrix, 
row  indexing, 
column  indexing, 
real  or*  imaginary  indexing, 
value  along  the  X-coordinate . 
value  along  the  Y-coordinate. 
dipole  number  indexing, 
radius  of  dipoles  . 

percentage  error  displacement  in  X(5). 
percentage  error  displacement  in  Y(5). 


OUTPUT  DATA: 

X(L) 

YCL) 

RAD 

PERRX 

PERRY 

CSSIM( IJ,KJ,N) 
IJ 
KJ 
N 

VLCI,K) 

I 

K 

TCOEFC I ,N , K) 


AMPLI(I,N) 

TCOEF 


X-coordinate  of  dipoles . 

Y-coordinate  of  dipoles . 
radius  of  dipole. 

percentage  error  displacement  in  X(5) 

percentage  error  displacement  in  Y(5) 

short-circuited  admittance  matrix. 

row  indexing. 

column  indexing. 

real  and  imaginary  indexing. 

phase  voltage  of  center  louded  dipole 

row  indexing. 

column  indexing. 

total  current  distribution,  where  I 
is  row  indexing,  N  is  column  indexing 
and  K  is  real  or  imaginary  indexing, 
amplitude  of  current  with  the  same 
indexing  as  TCOEF. 
is  also  punched  on  cards. 
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o  o  r>  o  non 


i  RCGR.AM  TC  CALCULATE  SHORT  CIRCUIT  ADMITTANCE 
PHASE  VOLTAGE  CF  CENTER  LOADED  DIPOLE 
AND  TOTAL  CURRENT  DISTRIBUTION  OF  ARRAY. 

DIMENSION  Y ( 5 , 5 ) , D I V ( 2 ) , YCCEF ( 5 , 2 ) »VL(5,2) 

DIMENSION  TRM15) »YTRM(5),P(1,2)  ' 

DIMENSION  ANG2{5) ,ANG3(5) ,ANG4(5),E(4,5,2) 

DIMENSION  CS(25,3G,2),CVI lO),COEF{10,10,2),CSSIM(lO,iO,2> 
DIMENSION  Cl (5,4,2) ,CCOEF{ 5,4,2) 

DIMENSION  TCI (30,4,2) ,TCOEF( 25, 4,2) 

DIMENSION  A V ( 2 ) 

DIMENSION  AMPL! (25,4) 

DIMENSION  X ( 5 ) , Y 1 ( 5 ) 

IPP  =  5 
NN  =  5 
t  H  = .  2  5 

AK=6. 2831853 
MM.=  I  PP+  l 
A*JN  =  NN 
DZ=EH/ANN 
I  N=2 
I0LT=5 

READ  IN  INVERTED  MATRIX 

READ (2, 101 ) ( ( (CS( I, J,K) ,K*i,2), J=l,30), 1=1,25) 

101  FORMAT ( 20A4  ) 

DO  102  1=1,25 
30  102  J  =  1 »  30 

in  102  k= l ,  2 

102  CS( I , J,K)=A4CVT(CS( I,J,K)  ) 

READ  IN  COORDINATES  CF  DIPOLES 
READ(2»250)(X(I),Y1(I)»I=1»NN) 

250  FORMAT! 2F0.4  ) 

C  PRINT  COORDINATES  CF  DIPOLES 

wR I TE (  5 , 254  > 

254  FORMAT ( 1H1 ,6X, ’COORDINATES  OF  DIPOLES*) 

CO  ?  5  3  I =1 ,NN 

255  WRI TE(5,25fi) I ,X( I ) , I ,Y1 ( I  ) 

2  56  FORMAT!  3X,  •X’,Il,'=*,F7.4,JX,*Y»,ti,,=',F7.4) 

C  READ  IN  RADIUS  OF  DIPOLE 

!<I  AC  (2,259)  RAP 

259  FORMAT ( F6. 4 ) 

WRITF  (:>,?60)RAU 

260  FORMA  T( 6 X,' RADIUS  OF  DIPOLE  IS',F7.4) 

KCA0(2,23T)PkRRX,  PERRY 

2  57  FORM AT ( 2F  7*4) 

WRIT£(5,258)PERRX, PERRY 

238  F0RMAK5X, 'PERCENTAGE  ERROR  DISPLACEMENT  IN  X5  AND  Y5  *  ,  /  15X,  •  X5=  • 
i,rb.4» 'Y5=* ,F8.4) 


C-2 


o  r>  o 


URIVING  COLUMN  MATRIX 

DO  3GO  M-ltMM 

AM=K 

Z={AM-1. )*C2 
RV=~SlNIAK*7)/60. 

300  CV ( X } =RV 

CALCULATE  SHORT  CKT  ADMITTANCES 

DO  701  KJ=i,NN 
JK=KJ*MM-MX 
DO  701  I J=1 » NN 
J=( I J-1)*IPP+1 
DU  702  K=l,2 
702  C0EF(J9KJsl  '=0.0 
DO  700  L=1,MM 
00  700  K=l,2 
LL= JK+L 

CHANG =CS  (  J  *  LI. » K  ) 

700  C0EF{ J,KJ,K)=CCEF{ J,KJ,K)+CHANG*CV(L) 

LO  701lK=i,2 

I  IK  =  -U 
N  =  1 1 K  +  3 

IF(N-1)20, 25,20 

25  CS$lM(lJ»KJ,N)=**COEF(J,KJ»IK) 

GO  TO  701 

20  CSSIM! I J,KJ,N)=COEF{ J,KJ, IK) 

701  CONTINUE 
hRITF.  (5,111) 

111  FORMAT! lHlt6X, ‘SHORT  CIRCUIT  ADMITTANCE  *  ) 

WR1TE( I  OUT, 921 ) ( ( (CSS  IX (1 J,KJtN) ,N= 1 , 2 ) , KJ= 1 , 5 ) , I J  =  l,5) 
921  FORMAT! 5X, El  3.6, 5X, El  3. 6) 

ANG 1=0.0 

RAD=3. 141593/130. 
f'EG  =  ANGl*RAf) 

C  CALCULATE  DRIVING  VCLTAGE  FOR  *0‘  PHASE 

1  =  1 

DO  100  K=l,/. 

C«)  105  J  =  l,4 
IF (K-l ) 3,3,5 
i  E( I , J,K)»CCS(CFG) 

GO  T(.  ,’u5 

5  L-  (  I  ,  J  ,K  )  =  S  Ii\  (  CEG  ) 

105  CUNT  I  NO! 

100  CONTINUE 

C  CALCULATE  DRIVING  VOLTAGE  FDR  ‘+90’  PHASE 

I  N  =  0 

DO  110  N=9C , 360 , 90 
F  N  =  I N  + 1 


C~3 


THETA=N-90 
OEG=THETA*RAD 
ANG2 ( IN ) =OEG 
110  CONTINUE 
1=2 

DO  Ub  K  =  1 » 2 
00  120  J  =  1 » A 
IF{K-1)9»7,9 
7  A=ANG2 ( J  } 

E ( I » J « K ) =COS ( A  • 

GO  TC  120 
9  8=ANG2 I J ) 

E{ I»J,K)=SIN(P) 

120  CONTINUE 
115  CONTINUE 

CALCULATE  DRIVING  VOLTAGE  FOR  »-90*  PHASE 
IN=0 

DO  125  N=9C,360,90 
IN=IN+1 
THE  T  A=90-N 
DEG=THETA*RAD 
ANG3 ( IN ) =OEG 
125  CONTINUE 
1  =  3 

DO  130  K=l,2 
DO  135  J  =  1  ,A 
IF(K-l)Ll»13»li 
13  A=ANG3 ( J ) 

E { I , J  »K ) =CCS { A ) 

GO  TO  135 
11  R  =  ANG3 ( J  ) 

E ( I »  J  t K ) =S IN ( 8 ) 

135  CONTINUE 
130  CONTINUE 

CALCULATE  DRIVING  VOLTAGE  FOR  *180‘  PHASE 
I  N  =  0 

DO  1 AO  NT=  1 »  2 
00  l A  5  N  =  1 80  *  360»  180 
I N= 1 N+  1 
THETA=N-180 
OEG=THETA«RAO 
ANGA (  IN ) =DEG 
IA5  CONTINUE 
1  AO  CONTINUE 
I  =A 

DO  150  K  =  1 » 2 
DO  155  J  =  1 »  A 
IF { K-l ) 15» 17  *  15 
W  A  =  ANGA ( J  ) 


El  I  *J*K)=COS( AJ 
GO  TO  155 
1 5  B=ANG4 l J ) 

E(  I » J,K)  =  SIMB) 

155  CONTINUE 
150  CONTINUE 

LOAD  ADMITTANCE 
YREAL=-G.Ol 
YIMAG=0.0 
I  J=5 

DO  160  KJ=1 , 5 
DO  160  K=l,2 

160  Y(KJ,K)=CSSIM{ IJ,KJ,K1 
CR£AL=Y(5,1) 

CIMAG=Y15,2) 

DIVli)=YREAL-CREAL 

DIV12)=-IYIMAG-C1MAG) 

D=DiV{l)**2+DlVi2)**2 
DO  170  1=1,4 
DO  180  K=1 ,2 
180  YCOEF ( I ,K ) =0.0 
00  175  J=l»4 
00  186  K=1 , 2 
TRMCKJ=E( I»J»K) 

185  YTRM(K)=Y{ J,K) 

CALL  CMULT ( TRM , YTRM, P ) 

DO  190  K=1 ,2 
YINCR=P l 1 , K ) 

190  YCOEF (I,K)=YCOEF(I»K) +YINCR 
175  CONTINUE 
170  CONTINUE 

DO  200  1=1,4 
DO  205  K=1 ,2 
TRM(K)=YCOEF( I ,K) 

205  YTRM(K)=DIV(K) 

CALL  CNULT ( TRM » YTRM , P ) 

DO  210  K= 1 , 2 
210  VL(I,K)=P(1,K)/D 

200  CONTINUE 

WR I TE ( 5, 1000 ) 

1000  FORMAT { 1H1 »6X, ’ PHASE  VOLTAGE  CF  CENTER  LOADED  DIPOLE* ) 
WR 1  T  E  ( 5 » 1 1 2 ) C (VLl I,KJ ,K=1,2) ,1=1,4) 

112  F0RMAT(2(5X,E13.6M 

J  =  5 

DO  333  1=1,4 
UO  333  K=1 ,2 
333  E(I,J,K>=VL(I,KI 
DO  5851=1,4 
DO  555J=1 , 5 


00  555  K=l,2 

555  CI(J,I,KJ=E(’»J,X) 
wRI7ciS,55?i 

557  FORMAT!///, 15X,*CI{J, I, K)  MATRIX*) 

WRITE (5, 556) H  iCl  I  J,ItK) ,K=1,2) , 1=1,4) ,J=1,S) 

556  F0RMAT(8(2X,FlG.i4) ) 

DO  320  1=1,5 

00  330  N=l,4 
DO  310  K=1 ,2 
310  CCCEF { I  ,N’,K }  =0.0 
DO  360  J=1 ,5 
DO  340  K=l,2 
TRH ( K )=CSSIM( I , J«K ) 

340  YTRM (K) =CI { J ,N»K) 

CALL  CMULT ( TRW, YTRH, P ) 

00  350  K=i , 2 

350  CCCEF ( I »N,k ) =CCOEF !I,N»K)+P(i,K) 

360  CONTINUE 
330  CONTINUE 
320  CONTINUE 

WR I TE  ( 5, 558  ) 

558  FORMAT!///, 15X , • CCGEF 1 1 ,N, K )  MATRIX  * ) 

WRITE (5, 556)  ( {  (CCOEFU  ,N,K )  ,K=1 ,2 )  ,N=1, 4 )  ,  1=1,5 ) 
C  CALCULATE  TOTAL  CURRENT  DISTRIBUTION 

C  CALCULATE  TCTAL  VOLTAGE  PER  PHASE 

DO  407  J  =  1 »  4 
DO  400  1=1,5 
DO  425  M=1,KM 
II»(I-1)*6+K 
AM=M 

Z=! AM-1. )*DZ 
AV { 1 ) =0 • 0 

AV ( 2 ) =~SIN{ AK*Z ) /60 . 

00  426  K=1 , 2 
TRM!K)=CI(I»J»K) 

426  YTRM { K ) =AV ! K  ) 

CALL  CMULT(TRM,YTRH,P) 

DO  427  K= 1 ,2 

427  TCI(II»J*K)=P(t,K) 

425  CONTINUE 

400  CONTINUE 
407  CONTINUE 

DO  450  1=1,25 
DO  450  N=l , 4 
DO  430  X=1 , 2 
430  TCCEF ( I,N,K)=0*0 
DO  450  J= ! , 30 
00  455  K=1 , 2 
TRM ( K ) =CS ( I , J , K ) 
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455  YTRM{K>=TCI( J,N,KJ 

CALL  CMlil  T { TRM , YTRH ,  P ) 

DO  450  K=J.,2 

450  TCQEF ( I  »N»K }  =TCGEF { I»N?K)-*-P(l»K) 

WRITE ( 5,666 ) 

666  FORMAT!///, 40X, ‘TOTAL  CURRENT  DISTRIBUTION* i 
DO  510  N=l,4 
DO  510  1=1,25 

510  AMPLIU,NJ  =  SQRT{(TC0EF(I,N,i))**2+(TC0EF(I,N,2)  )**2) 

*  WRIT£C5,99?*iUTC0£FU,N,K>,K=l,2),AHPLI  (I,N),N=l,4>,I*l,25) 

99?  FORMAT (6(1X,E13.6)) 

CALL  EXIT 


This  Fortran  II  program  plots  the  real  and  imaginary 
parts  of  the  current  vector  against  the  distance  along  the 
dipole.  Plotting  is  done  on  a  CALCOM  plotter. 


CUR(I,J,K)  - 

T 

j 

K 

X 


A  matrix  of  floating  point  numbers 

containing  the  array  to  be  plotted. 

row  indexing. 

column  indexing. 

real  and  imaginary  indexing. 

An  array  of  floating  point  numbers 
which  represents  the  distance  along 
the  dipole. 
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C  PROGRAM  TO  PLOT  CURRENT  DISTRIBUTION  OF  ARRAY. 

DIMENSION  w{6) ,X (6) 

DATA  X/O.OtC-. Co, 0.10, 0.1 5,0.20*0.25/ 

DIMENSION  C0R£Z5,A»2) 

READ  1 2,5) I ( I  CUR (I , J ,K ) , K=I ,2 ) , J= 1 ,4 ) , 1  =  1 ,25 ) 

5  FORMAT { AE13.6) 
rt{6)=0.C 
WRITE! 5,6) 

6  FORMAT  (1HD 
00  SOU  J=1 » A 

CALL  PS IZE ( 5.G , 10.0 ) 

CALL  PBOX 
CALL  PAXES 
DO  500  K=  l ,  2 
00  501  \  =  1,5 
DO  515  M  =  1 » 5 
I=(N-1>*5+M 
W{  M  )  =CUR  ( I  *  J  •»  K  ) 

515  WRITE (5*1)CLR{ I ,  J  ,  K  } ,  w  { M ) 

1  F0RVAT(2(2X,G13.6} ) 

CALL  PLCS  (0. 0, X, 0. 25,-0. 020, W, 0.020, 6) 

501  COMINUG 
PAUSF 

500  CON1 I NuF 
CALL  EXIT 
F.NC 
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This  Fortran  II  program  computes  the  input  to  output 
power  ratio  of  the  array  and  the  isolation  between  the 
center  element  and  the  circular  array. 

INPUT  DATA: 

A 
R 

X5 

Y5 
E 

V(I,J,K) 

CUR(I, J,K) 

OUTPUT  DATA: 

All  input  specified  above  are  written  out. 

PIN  -  input  power  to  each  dipole  on  the  circle. 

PO  -  output  power  to  the  fifth  element. 

PTOT  -  total  input  power. 

PTR  -  power  ratio, 

PDB  -  isolation  (in  DB). 


-  radius  of  dipoles. 

-  radius  of  the  circular  array. 

-  value  of  the  X-coordinate  of  the 
fifth  element. 

-  value  of  the  Y-coordinate  of  the 
fifth  element. 

-  percent  error  in  displacement  of  the 
fifth  element. 

-  voltage  sequence. 

-  current  distribution. 
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C  THIS  PROGRMA  COMPUTES  THE  INPUT  AND  OUTPUT  POWER  RATIO 
C  FOR  A  FIVE  ELEMENT  ARRAY, 

REAL  10 

DIMENSION  I0(2),CU(4,2) 

DIMENSION  CUR(25,4t2},V(4*4,2)f PIN(4, 2) ,P(4,2) 
DIMENSION  TRM(2) ,YTRM(2) ,REP(4,4) 

■C  IMPEDANCE 

ZL=100.0 
DO  150  1=1,4 
DO  150  K=l,2 
150  CUt l , K)=O.C 
C  REAO  IN  PHASE 

REA0(2,250)IFASE 
250  FORMAT (II) 

C  READ  IN  PARAMETERS  FOR  ARRAY 

READ(2,20)A,RTX5,Y5,E 
20  FORMAT  (4F6»4,F5«2) 

C  READ  IN  VOLTAGE  SEQUENCE 

READ t 2,230 ) ( ( ( V t I , J ,K ) , K=1 ,2 ) , 1=1 ,4 ) , J*l,4) 

230  FORMAT { 8F4 • 1 ) 

C  READ  IN  CURRENT  DISTRIBUTION 

READ (2,101) { { (CS( I,J,K),K=1,2), J=1 , 30) , 1=1 , 25 ) 

101  FORMAT { 20A4 ) 

DO  102  1=1,25 
DO  102  J= l , 30 
DO  102  K=1 » 2 

102  CS(  I,J,K)  =  A4CVT(CS( I , J , K ) ) 

DO  500  J=1 , IFASE 
PTOT=O.C 

DO  300  1=1,4 
1 1  =  (  I-l)*5  +  l 
DO  2b 1  K  =  1 »  2 

281  CUCI,KJ=CUR( II,J,K) 

C  CONJUGATE  CF  I 

CUR ( 1 1 »  J , 2 ) =-CUR ( II , J , 2 ) 

DU  325  K=l,2 
TRM(K)=CUR( 1 1 , J ,K ) 

325  YTRM(K)=V( I,J,X) 

CALL  CMULT(TRM,YTRM,P) 

DO  32?  N*l,2 
327  P IN ( I , N ) =P ( l ,N ) 

C  REAL  POWER 

REP ( I , J ) =P IN ( I , 1 ) 

PTCT=PTOT+REP( I, J) 

300  CONTINUE 

DO  305  N* 1,2 
305  1 0 ( N ) =CUR ( 2 1 , J ,N  J 

ABSI0*I0(1 )**2+I0(2)**2 
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P0=ABSi0*ZL 

PTR=PO/PTOT 

P0B=4.343*AL0G(PTR) 

WRITE(5,50) 

SO  FORMAT  (1HU 

WR I TE ( 5 » 355 ) J 

355  FORMAT ( 5X,  * PHASE=  *  f II*//) 

WRIT£(5,25)A,R»X5,Y5»E 

25  FORMAT (5X,'A=»,F8. 4, /, 5X, «R= • ,r 8.4 ,/,5X, *X5=' ,F8.4,/,5X, *Y5=' ,F8.4 
C,/,5X,'ERR0R=* ,F5.2,»  PERCENT' ,/// ) 
WRITE(5,400)({CU{I,K),K=l,2)«I=i,4) 

400  FORMATdX,  'Ii=',2F11.8,lX, »  I2=«  ,2FU.8»1X,«  I3=',2F11.8,1X,»I4=% 
C2FU.8,///) 

WRITE(5»405){ (V(I,J,K)tK»l,2JtI=l,4) 

405  FORMAT ( IX, • Vl  = ' , 2F11 .8 , IX, • V2= • , 2F11.8, IX, * V3= k  1 2F11.8, IX,  * V4=' , 
C2F11.8,///) 

wRITE(5,410) (REP( I,J) ,1=1,4) 

410  FORMAT (IX, • P 1= • , FI 1 . 7 , 2X , 'P2=«,F11.7,  2X, <P3=» ,FU.7,2XC »P4=« ,Fll. 
C7,//  ) 

*RITE(5,42Q)Zt,U0(I),I  =  l,2),ABSI0 
420  FORMAT ( IX, ' ZL  = ' ,  F6. 1 , 3X , ' 10= » , 2E11 .4 , 3X, • ABSI0=» , Ell .4, // ) 
WRITE(5»425)PT0T »P0 ,PTR»PD8 

425  FORMAT  (IX,  • PIN* * , Fli . 7,3X, ♦ P0= • , Fll . 7, 3X, • PO/P IN= • , F 11 .7, 3X , « PO/PI 
CN(D8)=' ,F11.7,///) 

500  CONTINUE 
CALL  EXIT 
END 
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This  Fortran  II  program  computes  and  plots  the  field 
pattern  and  phases  in  the  phi  plane. 

I 

INPUT  DATA: 


CURCI, J,K) 


I 

J 

K 

X(L) 

YCL) 


OUTPUT  DATA: 


A  matrix  of  floating-point  numbers 

containing  .the  current  at  the  res~ 

pective  nodes . 

row  indexing. 

column  indexing. 

real  or  complex  indexing. 

value  of  the  X-coordinate  of  the  L™ 

dipole . 

value  of  the  Y-coordinate  of  the  L™ 
cipole . 


IFASE 

CURCI, J,K) 
X(L) ,  YCL) 
FCOEFCNP,K) 
EMAGCNP) 

EDB 

PSHFTCNP) 


phase  sequence, 
as  above, 
as  above. 

complex  field  function, 
magnitude  of  FCOEF. 
value  of  normalized  EMAG  in  DB. 
phase  value  to  be  plotted. 
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C  PROGRAM  TO  CflMPl'Tfc  A  NO  PLOT  THE  FIELD  PATTERN  AND  PHASE 

C  OF  A  FIVE  ELEMCNT  CIRCULAR  ARRAY  IN  THE  PHI- PLANE 
DIMENSION  PhSFT  (  7 1>’) 

DIMENSION  SI (75) ,EMAG(75) ,X1 < 75) ,Y1 ( 75) 

Cl  MENS  I  (INCUR  (25,4,2)  , DCT(10),FASE(2)»  FCOEF  (75»2),X(5),Y(5) 
DIMENSION  TRM.(2  ) » YTRV ( 2 )  ,P { 1 ,2 ) 

REAL  Ml 
P I =3 . 14 1 59 
TWGPI=2.U*P1 
PI36=PI/36.0 
PI Z2=PI /2.C 
PI2=2.0*PI 
RAO=P  I / 180 .0 
DELTA=0.0 
EH= .25 
DZ  =  EH/5<>0 
THET  A=P  I Z2 
AS  I N  =  S!N { THETA ) 

ACCS=COS!THETA) 

C  READ  IN  CURRENT  AT  RESPECTIVE  NODES. 

READ (2,5) U (CUR(!,J,K)*K=l,2),J=i,A),I=l»25) 

4  F0RMAT(8(1X,E13.6) ) 

5  FORMAT ( 4E13.6) 

C  READ  IN  CORD INAT6S  OF  DIPOLES. 

REA0(2»10)(Xm»VU)»!»l*5) 

10  FORMAT { 2F6 . 4 ) 

C  READ  IN  PHASE 

DO  999  IFASE=1,4 
EMAX=0.0 

WR I TE ( 5 *  111 )  IFASE 
111  FORMAT (5X,tFASEs,I2) 

WRITE (5,4) {< (CUR( 1 ,J,K) ,K=1,2) ,J=1,4) ,1=1,25) 
WRITE(5,3)(X(I),Y(I), 1=1,5) 

3  F0RM.AT(2(IX,F8»4)  ) 

WRI TE ( 5 , 505 ) 

505  FORMAT ( 13X , ‘MAGNITUDE*  »15X» ’PHASE* ) 

DC)  2  1  =  1,73 
PHSFT ( I )=0.0 
on  2  J=1 ,2 
2  FCQFF ( I , J ) =0 .0 
DC  9  NP  =  i t  73 
PH  I  =  ( NP-1 ) *P I  36 

C  COMPONENTS  OF  UNIT  VECTOR  (A) 

AX=ASIN*COS(PHI ) 

AY=ASIN*SIN(PHI ) 

AZ= ACOS 
DO  24  1=1,5 
L  =  1 


015 


on  2‘j  LL 1=1 *9 
Ll=lO-LLl 
Ml  =  t>-Ll 

VECTOR  DOT  PRODUCT 

DOT  {  L )  =X  { I  )*AX+Y{  I)*AY+MiJS'DZ*AZ 

SOLID  ANGLE 

SI  { L ) =TkOP I *OOT ( L ) 

EASE ( 1 ) ?CGS l S I  {  L )  ) 
FAS6(2)=SIN{SI(L)} 

IF ( LL1-5 ) 6 » 14 » 7 


LL=-K1+ ( I-l ) #5+1 
GO  TU  11 
7  LL=Kl+(I-i)*5+l 

GO  TO  11 

14  LL=Ml+{ 1-1 )*5+l 
11  CONTINUE 

DO  30  N=1  ♦  2 

TRM(N)=CUR(LL, IF4SE,N) 

30  YTRM. ( N) =FASE (N ) 

CALL  CMULT ( TRM. » YTRM*  P ) 

DO  35  N= 1 , 2 

33  FCCEF(NP,N)=FCCEF(NP,N)+P(l,N) 

25  CONTINUE 
24  CONTINUE 

DO  36  N=l,2 

36  TRM(N}=FCaEF(NP,N) 

YTRM  i 1 ) =0 j  0 

YTRM ( 2 )=60.0*PI*S INI  THETA )*DZ 
CALL  CMULTI TRM,YTRM,P) 

00  37  N=1 f 2 

37  •  FCCEF(NP,N)=P(1,N) 

XI  (NP)  =  (NP-1)«PI36 
Z  =  FCQEF(NP,2)/FCf)EF(NP,l) 
IF(FCOEF(NPfl).GT.O)  GO  TO  38 
IF<FCOEF(NP,  :) .GT.O)  GC  TO  39 
1  =  3 

ARG=ABS( ATAN( Z ) ) 

GO  TO  50 
39  1=2 

ARG=PIZ2-A8S(ATAN(Z) ) 

GO  TC  50 

38  IF(FC0EF(NP»2) .GT.C)  GC  TO  41 

1=4 

ARG=PIZ2”A8S(ATAN(Z) ) 

GO  re  bJ 

41  1  =  1 

ARG=ARSCATAN(Z) ) 

30  PHSFT<NP)=AKG+U-n«PIZ2 

EMAG ( NP ) =SGRT ( FCOEF ( KP » 1 ) «*2  +  FC0EF ( NP, 2 ) **2 ) 


'  C-16 


V.R1  Te  ( 5, 510 )  EMAG (NP  )  ,PHL FT {NP  j 
510  FnRVAT{9X,E13.6,IlX,E13.fc) 

9  CONTINUE 

l*Rl  TE(2,5CC)  (EHAG(NP) , PHSPT  INP) , NP=1,72 ; 

500  FORMAT { 5X»8f 9.6) 

999  CONTINUE 

CALL  PSIZfciS. 0,5.0) 

CALL  PBGX 
CALL  PAXES 

CALL  PLOS  mG,Xl,TwCPI,Q.0,PHSFT,lQ.0,72) 
f:0  13  7  .VP=i  ,73 

EOC=2.68AO^ALOG { EMAG (NP ) /EMAX ) 

137  W<I TE( 5,902 ) (FCOEF(NP,I) ,1=1,2)  ,EMAG{ NP ) ,  ED6 
902  FORM AT(4 (5X, Eli. 6) ) 

CALL  PSIZE ( 5.G, 5.0 ) 

CALL  PBGX 
CALL  PAXES 
00  77  NP=1 , 7  ' 

XI  ( NP )  =  ( NP-  ,  , ‘  136 
Yl(NP)=EMAG(NP)/EMAX 
77  CONTINUE 

CALL  PLOTM G. 0, XI, TkCPI, 0.0, 71,1.0,72) 

PAUSE 
CALL  EXIT 
END 
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This  Fortran  II  program  computes  and  plots  the  field 
pattern  of  the  circular  array  in  the  Theta  plane.  (EMAG 
vs  THETA.)  Plotting  is  done  on  a  CALCOM  plotter. 

INPUT  DATA: 

CUR(I,J,K) 


I 

J 

K 

X(L) 

Y(L) 


OUTPUT  DATA: 

IFASE  -  phase  sequence. 

CUR( I, J ,K)  -  as  above. 

X(L),  Y(L)  -  as  above. 

FCOEF(NP,K)  -  complex  field  function. 

EKAG(NP)  -  magnitude  of  FCOEF . 

EDB  -  value  of  normalized  EMAG  in  DB , 


A  matrix  .of  floating-point  numbers 
containing  the  current  at  the  res¬ 
pective  nodes . 
row  indexing, 
column  indexing, 
real  or  complex  indexing, 
value  of  the  X-coordinate  of  the  L™ 
dipole. 

value  of  the  Y-coordinate  of  the  L"th 
cipole. 


C-18 


C  PROGRAM  to  compute  and  plot  the  field  pattern 
C  OF  A  FIVE  ELEMENT  CIRCULAR  ARRAY  IN  The  THETA  PLANE 
DIMENSION  SI  1755 ,EMAG ( 75) , XI { 75 ) ,Yi {75 ) 

DIM.ENSIOhCUR  ( 25*4 » 2 )  tOOT  (10)»FASE(2)  #  FCOEF  {75*2),X(5)»Y{5) 
DIMENSION  TRM(2) ,YTRM{2) ,P{ 1,2) 

REAL  Ml 

PI=3. 14159 

TW0PI=2.0*PI 

PI 36=P  1/36.0 

PIZ2=PI/2.0 

PI2=2.0«PI 

RA0=PI/180.0 

DELTA=OvO 

EH=.25 

DZ=EH/5.0  > 

PHI=0.0 
ACGS=C0S ( PHI ) 

ASIN  =  S IN ( PHI ) 

C  READ  IN  CURRENT  AT  RESPECTIVE  NODES. 

READ  (2*5)  (  (  (CUR(I»J,K),K--lt2),J=lt4),I  =  l,25) 

4  FORMAT(U{ IX.E13.6) ) 

5  FORMAT ( 4E l 3 . 6 ) 

C  READ  IN  CORCINATES  OF  DIPOLES. 

READ(2,lO)(X(I),Y(I)tI=l,5) 

10  FORMAT { 2F6.4 ) 

C  READ  IN  PHASE 

DO  999  IFASE=1»4 
EM4X=0.0 

WRITE(5»111 )  I  EASE 
ill  FORMAT { 5X» • FASE  =  '  12  ) 

WRITE { 5t4) { ( (CUR(ItJtK) , K=1 , 2 ) , J=1 ♦ 4 ) , 1= 1 ,25 ) 
WRITE(5f3)(X(I)fY(I)»Islt5) 

3  FORMAT (2{ IX , F8. 4 )  ) 

WR I  TE ( 5 » 901 ) 

901  FORMAT (9X, » REAL • , 15X ,  • I  MAG > , 1 5X ,  • MAG », 1 5X,  • DB » ) 

CO  2  1=1,37 
00  2  J  =  1 , 2 
2  FCOEF ( I , J ) =0 . 0 

DO  9  NP  =  i *  37 
THETA=(NP-1 

C  COMPONENTS  OF  UNIT  VECTOR  (A) 

AX=$IN{ THETA)*ACOS 
AY=SIN{ THETA)*ASIN 
AZ=COS{ THETA ) 

00  24  1=1,5 
1  =  1 

DO  25  LL1  =  1 »  9 
L I = 10-LL 1 
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-Vl=*>-L  l 

C  VECTOR  GOT  PRODUCT 

nCT(U=X{  IUA?X+Y{  l)$AY+Mi*DZ*AZ 
C  S«  it.  I D  AwGLE 

S I { L ) =T ,»UP  I  *UC  T  {  L  ) 

EASE { 1 ) =CGS ( SI ( L ) ) 
r ASE  C  2 ) =SIN ( SI tL) ) 

IF ( LL1-5 )6 , 14,7 

6  LL*-Fl+(  1-1  )*5  +  l 
GO  TO  11 

7  LL=.«1  +  U-1)*5-H 
GP  TO  11 

14  LL=*1+{I-1)*5+1 

11  CONTINUE 

DO  3C  N* 1.2 
TRf--{N)=CUR(LL,  IFASb,N) 

30  YTRM (N) =FASE (N ) 

CALI  CMoLT ( TRP , YTRM ,  P  ) 

PO  35  N  =  l,2 

3*3  FC0FF(NP»N)=FCCEF  ( NP ,N ) +P ( l » N ) 

25  CONTINUE 

24  CONTINUE 

00  3b  N=i,2 

,  3fc  TRNNI  =FCOEF  ( NP  »N  ) 

YTRM ( 1 ) =0. 0 

YTRM12)=60.C*PI*SIN(THETA)*U2 
CALL  CMULT{ TRP,YTRP,P) 

UO  37  N=l,2 

37  FCQ£F(NP,N)=P( l,N) 

EMAG{NP)=SCRT{FC0EF(NP,1)**2+FC0EF{NP,2)**2) 
IF  (EMAX-Ef'AG(NP) 1138,9,9 
138  EMAX=EMAG ( NP ) 

9'  CONTINUE 

00  137  N  P  =  1 , 3  7 

FDB=8.6860#AL0G(EMAG(NP)/ENAX) 

1 37  WRITE (5, 902) (FCOEFINP, I ) , I=i , 2 ) , EMAGI NP ) , ED8 
902  FnR"AT(4{5X,bl3.6) ) 

CALL  PSIZEtb. 0,5.0) 

CALL  PUOX 
CALL  PAXES 
00  77  NP=l »  37 
Xl(NP)=(NP-l)*PI3o 
Y 1 (NP)=EMAG{NP) /EKAX 
77  CONTINUE 

CALL  PLOTA(0.0,Xl,PI ,0.0, Y 1,1. 0,36) 

PAUSE 

999  CONTINUE 
CALL  EXIT 
END 
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APPENDIX  D 


THEORETICAL  DATA 

Data  that  were  n->t  included  in  the  text  are  appended 
here  for  reference.  They  include  plots  of  current  distri¬ 
butions  3  tables  of  short-circuited  admittances *  and  radia¬ 
tion  patterns . 


0.5%  POSITION  ERROR 
0°  PROGRESSION 


Figure  D-l:  Current  Distribution  for  the  C°  Phase  Progression 
Radius  CR)  of  the  Array  is  0.3 Radius  of  the 
Dipole  (A)  is  0.Q25X,  and  Position  Error  is  0.5%. 


RRENT 


0.5%  POSITION  ERROR 
180°  PROGRESSION 

- REAL 

-  IMAGINARY 

-  -  -  Jt  i 


\ 

N 


CURRENT 


Figure  D- 4:  Curren^  Distribution  for  the  0°  Phase  Progression. 

Radius  (R)  of  the  Array  is  0.3 ^  Radius  of  the 
Dipole  (A)  is  0.025  ^  and  Position  Error  is  10%. 
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Q%  POSITION  ERROR 
0%  -90°  PROGRESSION 

- -  REAL 

- —  IMAGINARY 

(  )  -90 5 


0,15  f  0,20  VO, 25 


Z(>0 


V(I2> 


for  the  +90°  and  -90°  Phase 
s  (R)  of  the  Array  is  0.3  A,  Radius 
0.025  X*  and  Position  Error  is  10%. 


CURRENT 


10%  POSITION  ERROR 
180°  PROGRESSION 

- REAL 

-  IMAGINARY 


FigureD-6:  Current  Distribution  for  the  180°  Phase  Progression 
Radius  (R)  of  the  Array  is  0.3^,  Radius  of  the 
Dipole  (A)  is  0.025  X,  and  Position  Error  is  10%. 


A  =  0.0250 
R  =  0.5000 
Xc=  0.0025 

y|-  o.oooo 

ERROR  =  0.5% 

Y11  =  0.12007E-01+j0.20060E-02 
Yl3  =  0 . 55746E-03+ jO . 91187E-04 
Yis  =  0 . 68221E-02+ jO . 37664E-02 
Y22  =  0 . 12103E~01+j0 . 19619E-02 
Y2u  =  0 . 55796£-G3+ jO . 90690E-04 
Y31  =  0 . 55746E-03+ jO .9118  6E-04 
Y33  =  0 . 12198E-01+j0 . 19161E-02 
Y35  =  0 . 70062E-02+j  0 . 36699E-02 
Y42  =  0 . 55796E~03+j0 . 90690E-04 
Y44  =  0 . 12Iu3E-01+j0 . 19619E-02 
Y51  =  0 . 6822lE-02+j0 . 37664E-02 
Y53  =  0 . 70062E-02+ j  0 . 36699E-02 
Y55  =  0 . 13993E-01+j0 . 11951E-03 

Zn  *  0 .11073E+03- jO  .  11507E+02 
Z13  =  0 . 26035E+02+j0 . 1957  GE+01 
Z15  =-0 . 6002 9E+Q2-jO . 859 10E+01 
Z22  =  0 . 11043E+03- j 0 .11714E+02 
Z24  =  0. 26033 E+02+j0, 19563E+01 
Z31  =  0 . 26035E+02+j0 .19670E+01 
Z33  =  0 , 11014E+03- jO , 11912E+02 
Z35  =-0 . 5967QE+02-j0 , 74185E+01 
Z42  =  0. 26033 E+02+jO, 19663E+01 
Z44  =  0.11043E+03-j0.11714E+02 
Z51  =-0.60029E+02~j0 , 85910E+01 
Z53  =-0.59  670E+02-j0.7418  5E+01 
Z55  =  0 . 14402E+03- jO .43584E+02 


Y12  =  0 . 3  3829E-02- jO . 3  5216E-03 
Y14  =  0 . 33829E-02- jO . 35216E-03 
Y2]  =  0 . 33829E-02-j0 . 3 5217 E- 03 
Y23  =  0. 34732E-02-j0.39713E-03 
Y25  =  0 . 69145E-02+ jO . 37184E-02 
Y32  =  0 . 34782E-02- jO , 3 9713 E- 03 
Y34  =  0 . 34782E-02-j0 . 39713E-Q3 
Y41  =  0 , 33829E-02- jO . 35217E-03 
Y43  =  0 . 34782E-02-jp  .  39713E-03 
Y45  =  0 . 69145E-Q2+jO ,  37184E-02 
Y52  =  0 . 69145E-Q2+ jO . 37184E-02 
Y54  =  0 . 69145E-G2+ jG . 37184E-Q2 


Z12  =-0 . 15186E+01+’)0 . 28999E+02 
Z14  =-0 . 15186E+01+j  0 , 28999E+02 
Z21  =-0 . 1518  6E+01+j0 . 28999 E+02 
Z23  =~0 . 18166E-f-014- j 0 .28796E+02 
Z25  =*-0 . 598  51E+Q2- jO  .  80004E+01 
Z32  s-0 . 18167E+Q1+ jO .  28796E+02 
Z34  =- 0.18 167 E+Ol+jO, 28796E+02 
Z41  =-0 ,15186E+01+j0 . 28999E+02 
Z43  =-0 .18166E+Gl+nQ . 28796E+02 
Z45  =-0 ,59851E+G2-jO .80004E+01 
Z52  =-0 , 59 8 51 E+02- j  0 . 80004E+01 
Z54  =-0.598  51E+02- jO .80004E+01 


Table  Dl:  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  (Z's)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-8 


A  =  0.0250 
R  =  G . 5000 
X5=  0.0100 
Y5=  0.0000 

ERROR  =2.0% 

Yn  =  0.117 2 5E-01+j0. 2117 5E-02  Y12 
Y 13  =  0. 55663E-03+j0.87410E-04  Yl4 
Y'lS  =  0. 65553E-02+j0 . 38881E-02  Y21 

Y  22  =  0 . 12109E-01+j0 . 19507E-02  Y23 
Y24  =  0. 56458E-03+j0.79495E-04  Y25 

Y  31  =  0. 556S3E-03+j0.87410E-04  Y32 
Y33  =  0.1248  5E-01+j0 . 17557E-02  Y34 
Y35  =  0 . 7  290SE-02+j0 .34996E--02  Y41 
Y42  =  0 , 56458E-03+jO . 79495E-04  Y43 
Y44  =  0. 12109E-01+j0.19507E-02  Y45 
Y51  =  0. 65553 E- 02 +j0. 38881E-02  Y52 
Y53  =  0 . 72905E-02+j0 .34996E-02  Y54 
Y55  =  0 . l4025E-01+j0 . 11909E-01 

Zn  =  0 . 11162E+03- jO  .  108  31E+02  Zl2 
Z13  =  0.2603  5E+02+j0 . 19800E+01  Z14 
Z15  =-0. 60528E+02-j0.10416E+02  Z21 

Z22  =  0,11042E+03-j0. 117 12 E+ 02  Z23 

Z24  =  0 . 26015E+02+j0 .19679 E+ 01  Z25 
Z31  =  0. 26035E+02+j0.19800E+01  Z32 
Z33  =  0 . 10924E+03- jO . 124  54E+02  Z34 
Z35  =-0 , 59095E+02- jO . 57250E+01  Z4i 
Z42  =  0. 26015E+02+j0.19680E+01  Z43 
Z44  =  0.11042E+G3-j0.11712E+02  Z45 
Z51  =-0. 60528 E+02-j0. 10416 E+ 02  Z52 
Z53  =-0 . 59095E+02- jO . 57249E+01  Z54 
Z55  =  0.14400E+03-j0 .43456E+02 


=  0.32439E-02-j0.29677E-03 
=  0 . 32439E-02- jO . 29677E-03 
=  0 . 32439E-02- jO . 29677E-03 
=  0. 36245E-02-j0.47782E-03 
=  0.69292E-02+n0 .36966E-02 
=  0 . 36245E-02- jO . 47782E-03 
=  0 . 3624 5E-02- jO . 47782 E- 03 
=  0. 32439 E-0 2- jO. 29677 E- 03 
=  0.36245E-02-jO .47782E-03 
=  0,69292E-02+j0.36966E-02 
=  0.69292E-02+j0.3  6966E-02 
=  0  .  ''9292E-02  + jO  .  36966E-02 


=-0.10830E+01+j0 . 29324 E+ 02 
=-0.10830E+01+j0.29324E+02 
=-0.10830E+01+j0.29324E+02 
=-Q . 22741E+01+ jO .  28  513E+02 
=-0.5983  5E+02- jO .80002 E+ 01 
=-0 .22741E+01+j0 .28513E+C2 
=-0 . 22741E+01+j0 . 28513 E+02 
=-0.10830E+01+j0 .29324E+02 
=-0.22741E+01+j0 .28513E+02 
=-0. 59835 E+02-j0. 80002 E+Oi 
=- 0,59835 E+0 2- jO .80002E+01 
=-0. 59835E+02- j0.80002E+01 


Table  D2:  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  (Z’s)  o^  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-9 


A  =  0.0250 
R  =  0.5000 
X5.-  0.0500 
Y5-  0.0000 

ERROR  =  10.0% 

Yll  =  0 . 3.033 6E-01+ jO .  23362E-02  Y12  =  0 . 25917E-02- jO .19167E-03 
Y13  =  0 . 52673E-03- jO. 62409E-05  Y14  =  0 . 25917E-02- jO .19167E-03 
Y15  s  0 . 53309E-02+ jO . 417  89E-02  Y21  =  0 . 25917E-02- jO . 19167E-03 
Y22  s  0 . 12253E-01+j  0 . 1S524E-02  Y23  =  0 .44076E-02- jO .12402E-02 
Y24  =  0 . 70792E-03- jO „ 2 18 8 IE- 03  Y25  =  0 . 72550E-02+ jO . 3 1139E-02 
Y31  =  0. 52673£-03-j0 . 62414E-05  Y32  *  0 . 440?6E-02-j0 .12402E-02 
Y33  =  0 . 13891E~01+j0 . 29074E-03  Y34  =  0 . 44076E-02- jO . 12402E-02 
Y35  =  0 . 88506E-02+j0 . 19367E-02  Y41  =  0. 25917E-02-j0 .19167E-03 
Y42  =  0. 70792E-03-j0 . 21881E-03  Y43  =  0 . 44076E-02- j 0 . 12402E-02 
Y44  =  0 . 12253E-01+ jO . 16524E-02  Y45  =  0 .72550E-02+j0 .31139E-02 
Y51  =  0. 53309E-02+j0.41789E~02  Y52  =  0 . 72550E-02+j0 . 31139E-02 
Y53  =  0 . 88506E-02+jQ . 19367E-02  Y54  =  0 . 72550E-02+ jO . 31139E-02 
Y55  =  0 . 14755E-01+ jO . 1Q760E-01 

Z11  =  0 .11583E+03-j0 . 57888E+01  Z3.2  =  0 . 94846E+00+j0 . 31453E+02 
Z13  =  0. 26046E+02+j0. 23150E+01  Z14  =  0 . 94846E+C0+j0 . 31453E+02 
Z15  =-0. 62163E+02-j0. 21522E+02  Z21  =  0 . 84846E+00+j0 . 31453E+02 
Z2o  =  0. 10995E+03-jO. 11666E+02  Z23  =-0 . 48784E+01+j0 . 27440E+02 
Z24  =  0. 25551E+02+j0, 20143E+01  Z05  =-0 . 59397E+02- jO . 79888E+01 
Z3I  =  0. 26046E+02+j0. 23150E+01  Z32  =-0 . 48784E+01+ jO . 27440E+02 
Z33  =  0. 10450E+03-j0. 14115E+02  Z34  a-0 . 48784E+01+j 0 . 27440E+02 
Z35  =~0 . 55428E+02+j0 . 20329E+01  Z41  =-0 .94S46E+00+ jO .31453E+02 
Z42  =  0. 25551E+02+jO . 20143E+01  Z43  =-0 .48784E+01+ jO .27440E+02 
Z44  =  0 . 10995E+03- jO  .  11666E+02  Zus  =-0 . 59397E+02-  jO  ,  79888E+01 
Z51  =-0.621o3E+02-j0. 21522E+02  Z52  =~0 . 59397E+02- jO  .79888E+01 
Z53  =-Q.55428E+02+j0 . 20330E+01  Z54  s-0 . 59397E+02- jO .7988 8E+01 
Z55  a  0«14338E+03-j0.40196E+02 


Table  D3 :  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  (Z's)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-10 


A  =  0.0250 
R  =  0.5000 
Xc=  0.0018 
Y5=  0.0018 

ERROR  =0.5% 

Yi:L  =  0 . 12Q34E-01+j0 . 19938E-02 
Y13  =  0 . 55772E-03+j0 .9Q923E-04 
Y t" c  =  0. 68481E-02+ jO  .  37  530E-02 
Y?2  =  0 .1203 4E-01+j0. 1993 3E- 02 
Y24  =  0 ,  55772E-03+ jO . 90923E-04 
Y31  =  0.5577  2E-03+j0 . 90924 E- 04 
Y33  =  0 . 12171E-01+ jO . 19290E-02 
Y35  =  0 . 69806E-02+i 0  .3  683  5E-02 
Y42  =  0 . 55772E-03+j0 .90924E-04 
Y44  =  0 .12171E~0I+j0 . 19290E-02 
Yc-,  =  0 . 68481E-02+j  0.37 530E-02 
Yc3  =  0 . 69806E-02+j0 . 36835 E- 02 
Y55  =  0.13993E-01+j0 .11951E-01 

Z1X  =  0 . 11065E+03- jO .11566E+02 
Z13  =  0. 26034E+02+jO .19666E+01 
Z15  =-0 . 59979E+02- jO .84248E+01 
Z22  =  0 . 11065E+03- j  0 . 11566E+02 
Z24  =  0 . 26034E+02+j0 . 19666E+01 
Z31  =  0 . 26034E+02+ jO .19666E+01 
Z33  =  0 . 11022E+03- jO .118  57E+02 
Z35  =-0 . 59721E+02- j  0 .7  5805E+01 
Z42  =  0 .26034E+02+j0 .19  66  6E+01 
Z44  =  0 .11022E+03- jO .11857E+02 
Z51  =-0.59979E+02-j0.84247E+01 
Z53  =-0.59721E+02~j0 .75804E+01 
Z55  =  0.14402E+03-j0.43584E+02 


Y12  =  0 . 33620E-02- jO .34266E-03 
Y14  =  0 . 34305E-02- jO .37432E-03 
Y2x  =  0 , 33620E-02- jO . 3426  5E-03 
Y23  =  0 . 34305E-02- jO .37432E-03 
Y25  =  0 . 68481E-02+ jC .37530E-02 
Y32  =  0 .34305E-02- jO .37432E-03 
Y34  =  0 .34992E-02- jO . 4 07 4 IE- 03 
Y41  =  0 .34305E-02- jO . 37432E-03 
Y43  =  0 .34S92E-02- jO .40741E-03 
Y45  =  0. 69806E-02+j0 .36835E-02 
Y52  =  0 . 68481E-02  +  n0 .37530E-02 
Y54  =  0.69506E-02+j0.36835E-02 


Z12  =~0-14533E+01+j0 .29045E+02 
Z14  =- 0. 16675E+01+j 0. 28897 E+ 02 
Z21  =~ 0.14 532E+01+j0. 29044 E+ 02 
Z23  =-0  ,1667  5E+01+ jO . 28897 E+ 02 
Z25  =-0.59979E+02-j0.84248E+01 
Z32  =-0 ,16675E+01+j0.28897E+02 
Z34  =-0 .18823E+01+ jO . 28753 E+02 
Z4{  =-0.1667  5E+01+ jO .28897 E+02 
Zu3  =-0.18823E+01+j0.28753E+02 
Z43  =-0.5872 IE ♦•02-jO  .7  5805E+01 
Z52  =-0 .59979E+02-j0 .84248E+01 
Z54  =-0 .59721E+02-j0 .75805E+01 


Table  D4:  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  (Z’s)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-ll 


0.0250 

0.5000 

0.0071 

0.0071 


A  = 
R  = 

X5  = 

V  - 


ERROR  =2% 

Yu  =  0 . 11837E-01+j0 . 2C721E-02 
Y13  =  0. 56063E-03+j0.83385E-04 
Y15  =  0 . 66652E-02+j0 . 38 3 3 IE- 02 
Y/'2  =  0 . 11837E-01+j0 . 20721E-02 
Y24  =  0 . 56063E-03+j0 . 83383E-04 
Y31  =  0. 56062E-03+j0.83385E-04 
Y33  =  0 . 12377E-01+ jO . 18149E-02 
Y35  =  0 . 7 187  IE- 02  + jO . 35570E-02 
Y42  =  0.56062E-03+j0.83382E-04 
Y44  =  0.123?7E-01+j0. 18149E-02 
Y51  =  0 . 6C652E-02+j0 . 38331E-02 
Y53  =  0 . 71871E-02+j0 .35570E-02 
Y55  =  0 . 14025E-01+j0 . 11908E-01 

Z11  =  0 . 11127E+0  3- jO . 11101E+02 
Z13  =  0 . 26025E+02+j0 . 19740E+01 
Z15  =-0 .60332E+02-j0.97010E+01 
Z22  =  0 . 11127E+03- jO . 11101E+02 
Z24  =  0. 26025E+02+j0 .19740E+01 
Z31  =  0 . 26025E+Q2+ jO  .  19740E+01 
Z33  =  0. 10 958 E+03-j0. 12252 E+02 
Z35  =-0.59313 E+02- jO. 63704E+QI 
Z42  =  0 .26025E+02+j0 .19740E+01 
Z44  =  0 .10958E+03- jO .12252E+02 
Z51  =-0.60332E+02-j0.97010E+01 
Z53  =-0 . 593 13 E+02- jO .63705E+01 
Z55  =  0. 14400 E+03-j0. 43455 E+02 


Yi2  =  0 . 3I649E-02- jO . 2 5 43 7 E- 03 
Y  X4  =  0 . 34334E-02- jO . 38186E-03 
Y21  =  0 .31649E-02- jO . 26436E-03 
Y23  =  0 . 34334E-02- jO . 38186E-03 
Y25  =  0 . 6S662E-02+j0 .38331E-02 
Y32  =  0 . 34334E-02- jO .38186E-03 
Y34  =  0 . 37052E-02- jO . 52156E-03 
Y41  =  0 . 34334E-02- jO .38186E-03 
Y43  =  0 .37052E-02- jO . 52156E-03 
Y45  =  0 .71871E-02+j0 .35570E-02 
Y',2  ~  0 ,66652E-02+j0 .38331E-02 
Y54  =  0 . 71871E-02  + jO . 3  5570E-02 


Z12  =-0.83570E+00+j0 .29509E+02 
Z14  =-0 , 167  61E+01+j0 . 28904E+02 
Z2i  =- 0.8 3569 E+OO+jO , 29509E+02 
Z-23  =-0.16761E+01+'j0  .28904E+02 
Z25  =-0 .60332E+02-j0.97010E+01 
Z32  =-0 . 16761E+01+ jO .28904E+02 
Z34  =- 0. 25266E+ 01 +j 0.28358 E+02 
Z41  =-0 . 16761E+01+ jO .28904E+0? 
Z43  =-0 . 25266E+01+j  0.28358 E+02 
Z45  =-0.59313E+02-j0.637Q5E+01 
Z52  =-0 ,60332E+02-j0 .97010E+01 
Z54  =-0 . 59313E+G2- jO .63704E+01 


Table  D5:  Short  Circuit  Admittance  (Y’s)  and  Open  Circuit 
Impedance  (Z*s)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-12 


A  =  0. 
R  =  0. 
X5=  0. 
V5=  0 

ERROR 

Y 11  = 
Y 13  = 
Yl5  = 
Y  22  = 
Y24  = 
Y31  = 
Y33  = 
Y35  = 
U2  = 

Y44  = 
YS1  = 

Y53  s 
Y55  = 

Zll  = 
Z 13  = 
Z15  =' 
Z22  = 
Z  24  = 
Z  3 1  = 
Z33  = 
Z35  = 
Z42  = 
Z44  = 

Zsi  = 
Z53  = 
Z55  = 


Table 


0250 
5000 
0354 
03  54 

=  10% 


0. 10912E-01vj0.22251E-02 
0 . 61713E-03- jO . 11477E-03 
0. 59290E-02+j0.38859E-02 
0. 10912E~01+j0.22251E-02 
0. 61713E-Q3-jO. 11477E-03 
0. 61713E-03-j0.11477E-03 
0.13443E-01+j0.73330E-03 
0. 84057E-02+j0.22736E-02 
0 . 61713E-03- jO . 11477^-03 
0.13443E-01+j0.73330^-0C 
0 . 59290E-02+ jO . 38859 E- 02 
0 . 84057E-C2+ jO .22736E-02 
0. 14742£-01+j0.10747E-01 


Y12  =  0. 22402 E-02-j0 . 1113 8 E- 03 
Y14  =  0.34899E-02-j0.58002E-03 
Y21  =  0 . 22402E-02- jO . Ill 38 E- 03 
Y23  =  0. 34899E-02-j0.58002E-03 
Y25  =  0 . 59290E-02+ jO . 38859 E- 02 
Y30  =  0. 34899 E-02-j0.58002E- 03 
Y34  =  0 . 47710E-02- jO . 16032 E- 02 
Y41  =  0. 34899 E-02-j0. 58002 E- 03 
Y43  =  0.47710E-02-j0.16032E-02 
Y4  £;  =  0 . 84G57E-02  + j  0 .22736E-02 
Y52  =  0.59290E-02+j0 . 388 59 E- 02 
Y54  =  0.84057E~02+j0. 22736 E- 02 


0 . 11418E+03- jO .795102+01 
0. 2580lE+02+j0.21559E+01 
0. 61582E+02-j0.17187E+02 
0 , 11418E+03- jO . 79511E+01 
0 .2580lE+02+j0 .21558E+01 
0. 2580lE+02+j0. 21558E+01 
0.1Q597E+03-j0, 13671 E+02 
0. 56646E+02-j0 .53209E+00 
0 ,2580lE+02+j0.21558E+01 
0 . 10597E+03- jO . 13671E+02 
0,61582  E+02- 'j  0.17 18  7E+02 
0 ,56646E+02-j0 .53204E+00 
0.14344E+03- j0,40199E+02 


7 12  =  0. 20759 E+ 01 +jO. 32659 E+02 
Z14  =-0.19004E+01+j0 .2S086E+02 
Z21  =  0. 20759 E+Ol+jO. 32659 E+02 
Z23  =-0.19003E+0I+j0.29086E+02 
Z25  -- 0. 61582 E+02-j0. 17187 E+02 

Z32  =~ 0.19 003 E+Ol+jO  .2908  6E+02 
Z34  =-0. 61323 E+01+j0,26939E+02 
Z41  =-0. 19 004E+01+j0. 29086 E+02 

Z40  -- 0. 61322E+ 01 +jO. 26939 E+02 
Z45  =-0.56846E+02-j0.53207E+00 
Z52  -- 0. 61582 E+02-j0. 17 187 E+02 
Z54  =-0. 56646E+02-j0.53209E+00 


D6 :  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  (Z*s)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 
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A  =  0.0250 
R  =  0.3000 
Xs=  0.0011 

y5=  0.0011 

ERROR  =0.5% 

',’n  =  0 . 61673E-Q2+j0 . 10064E-02 
Yi 3  =-0. 67796E-03-j0.42668E-03 
Yi5  =  0 . 74428E-04+ jO . 41473E-02 
Y22  =  0.6167  3E-02+ jO . 10064E-02 
You  =-G . 67796E-03- jO . 42668 E-03 
Y31  =-0. 677S6E-03-j0.42668E-03 
Y33  =  0.61738E-02+j0, 10484E-02 
v35  =  0. 72614E~04+j0 ,4 113 IE- 02 
Y42  =-0. 67796E-03- jO .42668E-03 
Y44  =  0 . 61738E-02+ jO . 10484E-02 
YS1  =  0.74430E-04+j0. 41473E-02 
Y53  =  0 . 72615E-04+ jO . 4 113 IE- 02 
Y55  =  0.45548E-03-j0.15182E-02 

Zn  =  0 . 12418E+03+j0 .91450E+01 
Z13  =- 0.153 43 E+02+j0, 3877  6E+02 
Z15  =  0 . 77707 E+ 01- j 0 . 53196E+02 
Z22  =  0. 12418E+03+j0 .91451E+01 
Z24  =-0. 15343E+02+j0. 38776 E+ 02 
Z31  =-0, 15343 E+02+j0.38778E+02 
Z33  =  0. 12457 E+03+j0 .90700E+01 
Z35  =  0.68860E+01-j0 .52572E+02 
Z42  =-0. 15343 E+02+j0. 3877 6E+02 
Z44  =  0. 12457 E+03+j0.90700E+01 
Z51  =  0.77707E+01-j0.53136E+02 
Z53  =  0. 69860 E+Ol-jO. 52572 E+ 02 
Z55  =  0.86780E+02+j0.55719E+02 


Y12  =  0.90189E-03+j0.12848E-02 
Y14  =  0 . 90513E-03+j0 .13059E-02 
Y21  =  0 . 90189E-03+ jO . 12848E-02 
Y2  3  =  0. 90513E-03  +  j0. 13059 E- 02 
Y2  5  =  0 .74426E-04  + jO .41473E-02 
Y32  =  0.80513E-03+j0 .13059E-02 

Y34  =  0 , 908  34E-03  + jO . 13268E-02 
Y41  =  0.90513E-03+j0.13059E-02 
Y43  =  0 . 90834E-03+ jO , 13268E-02 
Y45  =  0. 72613 E-04  +  j0 .4 113 IE- 02 
Y52  =  0. 74428 E-04+j0. 41473 E- 02 
Y54  =  0 . 7  2614E-04  + jO .4113  IE- 02 


Zi2  =-0. 50275E+02- jO  .33139E-*0: 
Zl4  =-0.50077E+02-j0.33175E+C: 
Z21  =-0.5027  5E+02- jO . 33139E+G2 
Z23  =-0.50077E+02-j0.33175E+0: 
Z  2  5  =  0 .77706E+01-j0 .53196E+02 
Z32  =-0. 50077 E+02-j0. 33175 E+0? 
Z34  =-0,498  81E+02- jO , 3  32 14  E+0 1 
Z41  =-0.5007r'E+02-j0  .33 17 5E +02 
Z43  =-0 .49881E+02-j0 .33214E+02 
Z45  =  0 , 698 59E+01- jO .52572E+02 
Z52  =  0 .77705E+01- jO . 53196 E+ 02 
Z54  =  0.69860E+01-j0.52572E+02 


Table  D7 :  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  (Z*s)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-14 


?=sra2S^UKB^«3SWiK5aK3«ir  ~**xr  X^.‘nmrr~r-rmv--*r-  *» 


M*as*v»aa*B>s£jaBC&&!aizs&S3fS£ asstes 


A  =  0.0250 
R  =  0.3000 
X5=  0.0042 
¥5=  0.0042 

ERROR  =2% 

Y 11  =  0.6158lE-02+j0.94689E-03 
Yi  3  =-0.  67789E-03- j0 ,  ‘‘2592E-03 
Yis  =  0 . 77093E~04+j0 . 41961E-02 
Y22  =  0. 6158lE-02+jG .94689E-03 
Y24  =-0. 67789E-03-j0.42592E-03 
Y  31  =-0.677  89E-03- jO . 42592E-03 
Y33  =  0 . 61827E-02+ jO , 11071E-02 
Y35  =  0 . 70126E-04+j0 . 40657E-02 
Y42  =-0. 67789E-03-j0.42592E-03 
Y44  =  0 . 61827E-02+ jO . 11072E-02 
Y51  =  0 . 7709lE-04+j0 , 41961E-G2 
Y53  =  0 . 7  0125E-04+ jO . 40657E-02 
Y55  =  0 . 45555E-03- jO .I5210E-02 

Zn  =  0. 123  62E+03+j0 .92371E+01 
Z13  =-0.15327E+02+j0 .38776E+02 
Z15  =  0 .88689E+01-jG.54045E+02 
Z22  =  0 . 12362E+03+j  0 .92371E+01 
Z24  =-0.15327E+02+i0.38776E+02 
Z31  =-0.15327E+02+j0 .38776E+02 
Z33  =  0 . 12512E+03+ jO .89501E+01 
Z35  =  0.58749E+01-j0.51665E+02 
Z42  =-0.15327E+02+j0.38776E+02 
Z44  =  0.12512E+03+j0.89502E+01 
Z51  =  0. 88689 E+01-j0,54045E+02 
Z53  =  0 . 53749E+01- j  0 . 5166  5E+02 
Z55  =  0. 86762 E+02+j0. 55634 E+02 


Yt  2  =  0 . 89263E~03+j0 . 12252E-02 
Yt4  =  0.90520E-03+30.13066E-02 
Y21  =  0.89263E-03+j0.122S2E-02 
Y  23  =  0 . 90520E-03+ jO .13066E-02 
Y25  =  0 .77091E-04+ jO .41961E-02 
Y32  =  0 .90520E-03+ jO .13Q66E-02 
Y34  =  0 .91728E-03+j0. 13855E-02 
Y41  =  0 . 90520E-03  +j 0  .13 066E-02 
Y43  =  0.91728E-03+jO. 13855E-02 
Y45  =  0 .70126E-04+j0 .40657E-02 
Y52  =  0 . 77091E-04+i0 .41961E-02 
Y54  =  0 . 70125E-04+j0 .40657E-02 


Z12  =-0.50835E+02-j0.33047E+02 
Z14  =-0 . 50061E+02- jO .33175E+02 
Z  21  =-0 ,50835E+02-j0 .33047E+02 
Z2  3  =-0 . 50061E+02-j0 .33175E+02 
Z2$  =  0 . 88688F+01- jO . 54 044 E+02 
Z32  =-0 . 50061E+02- jO .3317  5E+02 
Z34  =-0 . 49332E+02-jO .33333 E+02 
Z4X  =- 0 . 50061E+02- jO .3  317  5E+02 
Z43  =-0 . 49332E+02- jO . 3  3333E+02 
Z45  =  0.58748E+01-j0.51665E+02 
Z52  -  0.88688E+ 01- jO. 54044 E+02 
Z54  =  0.58748E+01-j0. 51665E+02 
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Table  D8  :  Short  Circuit  Admittance  (Y's)  and  Open  Circuit  .1 

Impedance  (Z’s)  of  Circular  Array  Elements  and 

the  Center  Element.  Values  Shown  Follow  the  1 

Order:  Real,  and  Imaginary.  ;■ 
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% 

1 


D^15 


A  =  0.0250 
R  =  0.3000 
X5=  0.0212 
Y5=  0.0212 

ERROR  =10% 

Yu  =  0 . 61027E~Q2+j0 . 61256E-03 
Yi3  =-0 . 67588E-03- jO . 40582E-Q3 
Y15  =  0 . 94486E-C4+j0 . 44801E-02 
Y22  =  0 . 61027E-02+jQ . 61256E-03 
Y24  =-0.67  587E-03- jO . 4Q582E-03 
Y31  =-0. 67587E-03-j0 .40582E-03 
Y33  =  0. 62296E-02+j0.14206E-02 
Y35  =  0. 57071E-G4+j0.3818GE-02 
Y42  =-0 . 67587E-03- jO . 40582E-03 
Y44  =  0 . 62296E-02+j0 . 14206E-02 
Y51  =  0.9i'486E-04+j0.44801E-02 
Y53  =  0 . 5707lE-34+j0 . 38180E-02 
Y55  =  0 ,45785E-03-j0 . 15949E-02 

Zn  =  0.12052E+03+j0. 94596E+01 
Zi 3  =-0 . 14910E+02+j0 . 38778E+02 
Z15  =  0.14o39E+02-j0. 58051E+02 
Z22  =  0 . 12052E+03+j0 . 94597E+01 
Z24  =-0.14910E+02+j0 .38778E+02 
Z31  =~0.14910E+02+j0 .38778E+02 
Z33  =  0,12 794 E+03*j0 . 80123E+01 
Z35  =-0 . 22818E+00- jO .46143E+02 
'42  =-0. 149 10E+02+j 0.38778 E+02 
Z44  =  0 . 12794E+03+j0 .80123E+01 
Z51  =  0 . 14 63 9 E+02- jO .58051E+02 
Z53  =-0 , 22821E+00-j0 . 4 614 3 E+02 
Z55  =  0 , 86315E+02+jQ .53452E+02 


Y12  =  0 .83726E-03+j0 . 89092E-03 
Y14  =  0 . 90721E-Q3+j0 .13267E-02 
Y21  =  0.83727E-03+j0.89092E-03 
Y23  =  0 .90722E-03+jO .13267E-02 
Y25  =  0 . 94485E-04+j0 . 44801E-02 
Y32  =  0 . 90722E-G3+j0 . 13 2 67 E- 02 
Y34  =  0 .9  6417E-03+j0 . 16990E-02 
Y41  =  0 ,90722E-03+j0 .13267E-02 
Y43  =  0 ,96417E-0?+j0 .16990E-02 
Y45  =  0 ,57071E-04+j0 .38180E-02 
Y52  =  0 . 94485E-04+ jO . 44801E-G2 
Y54  =  0 . 57071E-04+j0 . 38180E-02 


Z12  =-0 . 53936E+02-j0 .32824E+02 
Z14  =-0 .49 64 3 E+02- jO .33173E+02 
221  =-0.53936E+02-j0.s2824E+02 
Z23  =-0  ,49643E+02-j  0 . 33173E+02 
Z25  =  0 . 14639E+02- jO . 58051E+02 
Z32  =-0 .49 64 3 E+02- jO .3317 3 E+02 
Z34  =~0 .46506E+02- jO . 34271E+02 
Z41  =-0 ,49643E+02-j0 .33173E+02 
Z43  =-0.46506E+02-jQ.34271E+02 
Z45  =-0 . 22827E+00- jO .46143E+02 
Z52  =  0 ,14639E+02-j0 . 58051E+02 
Z54  =-0 . 22827E+00-j0 .46143E+02 


Table  D9:  Short  Circuit  Admittance  (Y’s)  and  Open  Circuit 
Impedance  (Z's)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-16 


A  =  0.0067 
R  =  0.3000 
X5=  0.0060 
Y5=  0.0000 

ERROR  =  2% 

Yjt  =  0 . 65745E-G2- jO . 29000E-02  Y12  =  0 ,97910E~03+j0 .13669E-02 
v-,3  =-0 , 72367E-03- jO  . 40158K-03  Y1U  =  0 .97910E-03+j0 . 13669E-Q2 
Yls  =  0 . 14269E-03+j0 . 45175E-02  Y2i  =  0 .97910E-03+j0 .13669E-02 
Y22  =  0.65960E-G2-j0.27761E-02  Y93  =  0 . 10003E-02+ jO .14868E-02 
Y2u  =-0 . 72399E-03-iQ . 40029E-G3  Y25  =  0 . 13623E-03+j0 .44180E-02 
Y31  =-0 . 72367E-03- jO . 40158E-03  Y32  =  0 .10003E-02+ jO . 14868E-02 
Y33  =  0 . 66170E-02- jO . 26600E-02  Y34  =  0 . 10003E-02+j0 . 14868E-02 
Y35  =  0 . 12928E-03+j0 . 43226E-02  Y41  =  0 .97910E-03+j0.13669E-02 
Y42  =**0 . 72399E-03~j0 . 40029E-03  Y43  =  0 .10003E-02  + jO  .14858E-02 
Y44  =  0 . 65960E-02-j0 . 27761E-02  Y45  =  0 . 13623E-03+ jO . 44180E-02 
Y51  =  0 . 14269E-03+j0 . 45175E-02  Y52  =  0 . 13623E-G3+ jO .44180E-G2 
Y5J  =  0. 12928 E-03+i0  .IJ3226E-02  Y54  =  0 . 13623E-03 + jO  .44180E-02 
Y55  =  0 . 48506E~03-j0 . 53367 E-02 

Zn  =  0 . 8  8838E+02+j0 . 37  370E+02  Z12  =-0 . 1C606E+02-  jO  .41101E+02 
Zi 3  =-0 . 34531E+02- jG . 2S428E+01  Z14  =-0 .10606E+02- jO .41101E+02 
Z15  =  0 . 24719E+02- jO . 41721E+02  Z21  =-0 . 10606E+02-j0 . 41101E+02 
Z22  =  0,89063E+02+'j0.37573E+02  Z29  =-0 .10379E+02-j0  .40911E+02 
Z<>4  =-0 . 34530E+02- jO .  25418E+01  Z25  =  0 .22761E+02- jO  .  41855E+02 
Z31  =~0 . 34531E+02~j  0.25428E+01  Z39  =-0 .10379E+02-j0 .40911E+02 
Z33  =  0 . 8929lE+02+j0 . 37750 E+ 02  Z34  =-0 . 10379E+02-i0 . 40911E+02 
Z35  =  0 . 20856E+02~jO .41918E+02  Z41  =-0 . 10606E+02- jO . 4I101E+02 
Z49  =-0.34530E+02-j0.25418E+01  Z43  =-0 .10379E+02-j0 .40911E+02 
Z44  =  0.87063E+02+j0. 37573E+02  Z45  =  0 . 227 61E+02- j 0 .41855E+02 
Z51  =  0.24719E+02-j0.41721E+02  Z52  =  0 . 22761E+02- jO . 41855E+02 
Z53  =  0.208  56E+02- jO .41918 E+ 02  Z54  =  0 ,227 61E+02- jO  .41855E+02 
Z55  =  0. 74842 E+ 02 +j0 . 39682 E+ 02 


Table  DIO:  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 
Impedance  CZ’s)  of  Circular  Array  Elements  and 
the  Center  Element.  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-17 


A  =  0.0123 
R  =  0.3000 
Xr=  0.0060 
Y5=  0.0000 

EFROR  =2% 

Yu  =  0 . 6398  0E- 02- jO  .  19  50  IE- 02 
Y 13  =-Q . 70168E-03-i0 . 39097E-03 
Y15  =  0 .14172E-03+j0  .43936E-02 
Y  22  =  0 . 64 19 IE- 02- jO .18293 E-02 
Y24  =-0 . 701S9E-03-j0 . 38970E-03 
Y31  =-0.70268E-03-j0.39097E-03 
Y33  =  0 . 64397E-02- j  0 . 17161E-02 
Y35  =  0. 12908E-03+j0 . 42038E-02 
Y4?  =-0 . 70199E-03-j 0 .38970E-03 
Y44  =  0 . 64191E-02- jO . 18293E-02 
Y51  -  0 . 1417  2E~03+j0 . 43936E-02 
Y53  =  0 . 12908E-03+j0 .42 038 E-02 
Yss  =  0 . 47422E-03-j0.43123E-02 

Z11  =  0.97609E+02+j0.34704E+02 
Zl3  =-0 . 36886E+02+j0 . 76726 E+01 
Z15  =  0. 21457 E+02-j0 . 45055E+02 
Z22  =  0.98030E+02+j0.34948E+02 
Z24  =-0.36884E+02+j0.76739E+01 
Z31  =-0 . 3  6386E+02+ jG .76726E+01 
Z33  =  0.98445E+02+j0 .35147E+02 
Z35  =  0.17190E+02-j0.44560E+02 
Z42  =-0.36884E+Q2+j0.76739E+01 
Z44  =  0.98030E+02+j0.34948E+02 
Z51  =  0. 21457 E+0 2- jO. 45055 E+02 
Z53  =  0 . 17190E+02- jO .44560 E+02 
Z55  =  0 . 76077E+02+j0 .42386 E+02 


Y12  =  0.95574E-03+j0 .13 2 28 E-02 
Y14  =  0 .95574E-03+j0  .3.3228E-02 
Y21  =  0 .95574E-03+j-0 .13228E-Q2 
Y23  =  Q.97660E-03+j0.i4398E-02 
Y25  =  0 .13565E-03+ jO .42967E-02 
Y32  =  0 .97661E-03+ jG .14398E-02 
Y34  =  0 .97661E-03+ jO .14398E-02 
Y41  =  0.95574E-03+j0.13228E~02 
Y43  =  0.97660E-03+j0.14398E-02 
Y45  =  0. 13 565E-03+j0. 42967 E-02 
Y52  =  0 .13565E-03+j0 .42967E-C2 
Y54  =  0. 13 565 E-03+j0. 42967 E-02 


Z12  =-0  .17955E+02-j  0.42433 E+02 
Zl4  =-0 .17955E+02-j0 .42433E+02 
Z21  =-0,17 955 E+0 2- jO. 42433 E+02 
Z23  =-0 ,17537E+n2-j0.42212E+02 
Z25  =  0.19 286 E+02-j0 .44846E+02 
Z32  =-0 .1753 7E+C 2- jO .42212 E+02 
Z34  =-0 .17537E+02-j0.42212E+02 
Z41  =-0 .17955E+02-i0.42433E+02 
Z43  =- 0.17 537E+0 2- jO. 42212 E+02 
Z45  =  0.19286E+02-j0 .44846E+02 
Z52  =  0 .19286E+02-j0 .44846E+02 
Z54  =  0.19286E+02-j0.44846E+02 
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Table  Dll:  Short  Circuit  Admittance  (Y!s)  and  Open  Circuit 
Impedance  (Z*s)  of  Circular  Array  Elements  and 
the  Center  Element*  Values  Shown  Follow  the 
Order:  Real,  and  Imaginary. 


D-18 


''4 
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1 


A  =  0.0135 
R  =  0.3000 
X5=  0.0080 
Y$  =  0.0000 

ERROR  =2% 

Yi  i  =  0 . 35926E-03+ jO . 10541E-01 
=-0.13368E-03-j0.45604E-04 
Y15  =  0 . 12006E-03- jO .19757E-03 
Y22  =  0. 35988E-03+j0.10541E~01 
Yon  =-0.13367E-03-j0.45608E-04 
Y3l  =-0 . 13368E-03- jO .45604E-04 
Y33  -  0 .3o049E-03+j0 . 10542E-01 
Y3r,  =  0.10577E-03-j0. 19656E-03 
Y42  =-0 . 13367E-03- jO . 45608E-Q4 
Y44  =  0.35988E-03+j0. 10541E-01 
Y51  =  0 . 12006E-03- jO . 19757E-03 
Y53  =  0 . 10577E-03- jO . 19656E-03 
Y55  =  0 . 32123E-03+j0 . 10546E-01 

Zn  =  0. 32539 E+Ol-jO, 94822E+02 
Z13  =-0.11370E+01-j0. 58803E+00 
Z.15  =  0 .12124E+01-j0.17  881E+01 
Z22  =  0 . 32569E+01- j»; .  94819E+02 
Z04  =-0 . 11370E+01- jO . 56805E+00 
Z?i  =-0.11370E+01-j0.56803E+00 
?33  =  0. 32599 E+Ol-jO. 94817 E+ 02 
Z35  =  0 .1084UE+01- jO .17887E+01 
Z49  =-0.11370E+01-j0.56805E+00 
Z44  =  0. 32569 E+Ol-jO. 94819 E+ 02 
Zt;->  =  0.1212“E+Ci-j0.17881E+01 
Z53  =  0.10844E+01-j0.17887E+01 
Z55  =  0. 30503 E+ 01- jO. 94813 E+ 02 


Y12  =-0.20373E-04-j0.18020E-03 
Y14  =-0.20373F-04-j0 .18020E-03 
Y21  =-0 . 20373E-04- jO .18020E-03 
Y23  =~0 .19759E-04- jO . 17983E-03 
Y25  =  0 .11284E-03-j0 . 19710E-03 
Y32  =~0 .19759E-04-j0 . 17953E-03 
Y34  =-0 .197  59E-04-j0  .17983E-03 
Y41  =-0 . 20373E-04- jO . 18 020 E- 03 
Y43.  =-0 . 19759E-04- jO  .17  983E-03 
Y45  =  0 . 11284E-03- jO .19710E-03 
Y52  =  0 . 11284E-03- jO . 19710E-03 
Y54  =  0.11284E-03-j0.19710E-03 


Z12  =-0 ,7088SE"0x- j 0 .16640E+01 
Z14  =-0  ,70889E-01-j0 . 16640E+01 
Z21  =-0 .70888E-01-j0 .16640E+01 
Z23  =-0.67864E-01-j0.16G17E+01 
Z2S  =  0 .11478E+01- jO . 17887E+C1 
Z32  =-0 ,67862E-01-j0.16617E+01 
Z34  =-0 . 67862E-01- jO . 16617 E+ 01 
Z41  =- 0. 70888 E- 01- jO. 16640 E+ 01 
Z4 *3  =-0 .67864E-0i-j0 .16617E+01 
Z45  =  0 .11478E+01-j0 .17887E+01 
Z52  =  0 . 11477E+01-i0 . 178  87E+01 
Z54  =  0 . 11478E+01- jO .178  87E+01 


abK  T12:  Short  Circuit  Admittance  (Y's)  and  Open  Circuit 

Impedance  (Z's)  of  Circular  Array  Elements  and  the 
'"‘enter  Element.  Values  Shown  Follow  the  Order: 
Peal,  and  Imaginary. 
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FR  (db) 


RfLATIVF.  POv,hH  (db! 


A  =  0,0123  X  R-  0,3000  X 
X5=  0,0060  X  Yj=  0,0000  \ 
ERROR  =  2,0% 


9 


i 


$  -  PUNE  ANGLE  (Decrees) 

figure  D-8  :  Theoretical  Azimuthal  Plane  Radiation  Pattern 
of  the  Circular  Array. 
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RELATIVE  .\flO  (ee) 


A  =  0,0067  A  R  =  0,3000  A 
Xf  0,0060  A  Ys=  0,0000  A 
ERROR  =  2.0% 


n 

* 

-1 

-6 

-10 


$  -  PLANF  ANGLE  (Degrees) 

Figure  D-1Q  :  Theoretical  Azimuthal  Plane  Radiation  Pattern 
of  the  Circular  Array. 
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